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PKEFACE. 



The object of the following treatise is to exhibit, in a concise 
form, the elementary properties of the expressions known by 
the name of Laplace's Functions, or Spherical Harmonics, 
and of some other expressions of a similar nature. I do not, 
of course, profess to have produced a complete treatise on 
these functions, but merely to have given such an introduc- 
tory sketch as may facilitate the study of the numerous 
works and memoirs in which they are employed. As 
Spherical Harmonics derive their chief interest and utility 
from their physical applications, I have endeavoured from 
the outset to keep these applications in view. 

I must express my acknowledgments to the Rev. C. H. 
Prior, Fellow of Pembroke College, for his kind revision of 
the proof-sheets as they passed through the press. 



N. M. FERRERS. 



GONTILLE AND GaiuS CoLLSOS, 

August, 1877. 
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INTRODUCTORY. DEFINITION OP SPHERICAL HARMONICS. 

1. If V be the potential of an attracting mass, at any 
point X, y, z, not forming a part of the mass itself, it is 
known that Fmust satisfy the differential equation 

d^V d^V d^V ^ 

^+^+^=« «' 

or, as we shall write it for shortness, V*F=0. 

The general solution of this equation cannot be obtained 
in finite terms. We can, however, determine an expression 
which we shall call F,, an homogeneous function oix,t/, z 
of the degree % % being any positive integer, which will 
satisfy the equation ; and we may prove that to every such 
solution V^ there corresponds another, of the degree — (i + 1), 

expressed by -^ , where r* = a?" + y^ + «*. 

For the equation (1) when transformed to polar co-ordi- 
nates by writing a; = r sin ^ cos ^, y = r sin ^ sin 0, « = r cos ^, 
becomes ' 

And since F satisfies this equation, and is an homo- 
geneous function of the degree ij F, must satisfy the equa- 
tion 

t (i + 1) F, + -v— . -771 sm ^ -^* + -^.-A-^ = 0, 
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since this is the form which equation (2) assumes when V 
is an homogeneous function of the degree u 

Now, put T^ = r^*** Utf and this becomes 
or 

Now, since t^ is a homogeneous function of the degree 

-(i + 1), 

dr * dr 

dr* dr 

therefore equation (2) becomes 



d^irV) ^ 1 d (.M\ , 1 cZ^t^_ 



air siti ^ c/^ 



shewing that U^ is an admissible value of F, as satisfying 
equation (2). 

It appears therefore that every form of U^ can be ob- 
tained from V^y by dividing by r'^^y and conversely, that 
every form of F< can be obtained from U^ by multiplying: 
by r***\ Such an expression as FJ we shall call a Solid 
Spherical Harmonic of the degree i. The result obtained 
by dividing F< by r\ which will be a function of two inde- 
pendent variables ^ and <^ only, we shall call a Surface 
Spherical Harmonic of the same degree. A very important 
form of spherical harmonics is tViat \<Yv\c\i as mde^endent 
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of ^. The solid harmonics of this form will involve two of 
the variables, x and y, only in the form a?" + y", or will be 
* functions of a?* + y* and z. Harmonics independent of (j) are 
called Zonal Harmonics, and are distinguished, like spherical 
hamonics in general, into Solid and Surface Harmonics. 
The investigation of their properties will be the subject of 
the following chapter. 

The name of Spherical Harmonics was first 'applied to 
these functions by Sir W. Thomson and Professor Tait, in 
their Treatise on Natural Philosophy. The name " Laplace's 
Coefficients" was employed by Whewell, on account of 
Laplace having discussed their properties, and employed 
them largely in the Mecanique Celeste. Pratt, in his 
Treatise on the Figure of the Earth, limits the name of 
Laplace's Coefficients to Zonal Harmonics, and designates 
all other spherical harmonics by the name of Laplace's 
^m£on8. The Zonal Harmonic in the case which we shall 
consider in the, following chapter, i.e., in which the system 
is symmetrical about the line from which 6 is measured, 
^as really, however, first introduced by Legendre, although 
t^e properties of spherical harmonics in general were first 
<ii8cu8sed by Laplace; and Mr Todhunter, in his Treatise, 
<^n this account calls them by the name of " Legendre's 
Coefficients," applying the name of "Laplace's CoeflScients" 
to the form which the Zonal Harmonic assumes when in 
place of cos^, we write cos ^cos^' + sin dsin^'cos(<^ — <^'). 
?^fae name " Kugelfunctionen ** is employed by Heine, 
iQ Ms standard treatise on these functions, to designate 
Spherical Harmonics in general. 



\— ^ 



CHAPTER IL 



ZONAL HARMONICS. 



1. We shall in this chapter regard a Zonal Solid Har- 
monic, of the degree i, as a homogeneous function of 
{pi? +y*)i, and z, of the degree % which satisfies the equation 

dji?'^ d/'^ dz^'^ ' 

Now, if this be transformed to polar co-ordinates, by 
writing r sin 6 cos <\> for a?, r sin d sin ^ for y, r cos for ;?, we 
obser\'e, in the first place, that a^+y* = r*sin*^. Hence 
V will be independent of ^, or will be a function of r 
and ctoly. The differential equation between r and 
which it must therefore satisfy will be 

d?{rV) ^ 1 d f . .dF> 



<^7^ 






Now y, being a function of r of the degree % may be 
expressed in the form r^P^y where P^ is a function of only. 
Hence this equation becomes 

or, putting cos O^fi, 

In accordance with our definition of spherical surface 
harmonica, P^ will be the zonal surface '^aaxTSioma of the 
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degree i. When it is necessary to particularise the variable 
involved in it, we shall write it P, (jii). 

The line from which is measured, or in other words 
for which /Lt = l, is called the Axis of the system of Zonal 
Harmonics; and the point in which the positive direction 
of the axis meets a sphere whose centre is the origin of 
co-ordinates, and radius unity, is called the Pole of the 
system. 

Any constant multiple of a zonal harmonic (solid or 
surface) is itself a zonal harmonic of the same order. 

2. The zonal harmonic of the degree i, of which the 
line /Lt = 1 is the axis, is a perfectly determinate function of 
/i, having nothing arbitrary but this constant. For the 
expression r*F^ may be expressed as a rational integral 
homogeneous function of r and ^, and therefore P, will be 
a rational integral function of cos 0, that is of /i, of the 
degree i, and will involve none but positive integral powers 

of /A. 

But P, is a particular integral of the equation 



i {('-"•) 



•*./(")) 



clfJb 



J 



+i(i-+i)/M=o.:....(3), 



and the most general form of f(jj) must involve two ar- 
bitrary constants. {Suppose then that the most general 

form o{f{/M) is represented by P, Ivdfi. We then have 

+2a-/'-)f.H-i'.|{(i-M'i..}. 

Hence, adding these two equations together, and ob- 
serving that, since P^ satisfies the equation ^^S^,\!ti^ e^^^'SkJsvsec^. 
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of jvdfi will be identically eq[ual to 0, we obtain, for the de- 
tennination of v, the equation 

wiience .P,(l-M*)g+2|(1-/*') ^-/*P.)»^0, 



or 



?^(if'-i^-)^-»- 



the integral of which is 

log » + log P* (1 - fi*) = log (7j = a constant ; 

_ c^, 

Hence Jt)d/i ^C+C.j p» n -u'j 5 

and we obtain, for the most general form off(jj,), 

* Now, P^ being a rational integral function of /t of i 

dimensions, it may be seen that i,^ ^rrm will assume the 

form of the sum of i + 2 logarithms and i fractions, and 
therefore cannot be expressed as a rational integral function 

of fjb. Expressions of the form P^ j j- ^rwi ^^e called KugeU 

functionen der zweiter Art by Heine, who has investigated 
their properties at great length. They have, as will hereafter 
be seen, interesting applications to the attraction of a sphe- 
roid on an external pomt. We shall discuss their properties 
more fully hereafter. 

. 3. We have thus shewn that the most general solution 
of equation (2) of the form of a rational integral function of a 
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involves but one arbitrary constant, and that as a factor. 
We shall henceforth denote by P„ or P, (/i), that particular 
form of the integral which assumes the value unity when /a 
is put equal to unity. 

We shall next prove the following important proposition. 

l/hbe less than unity y and if (1 - 2/ih + h*)"^ he eocpanded 
in a series proceeding hy ascending 'powers of h, the coeffioient 
ofljiwiUheV^. 

Or, (l-2/Ah + hT* = P,+ P> + ... + P,h*+... 

We shall prove this by shewing that, if H be written for 
(1 — 2/xA + A*)"^, H will satisfy the dififerential equation 

For, since H=(l-2iih + hy\ 

I dff _ jn ^ 

h d/i ' 

= - 2/*ff' + 3(1- fi*) H' '^J- 
= -2/*a» + 3(l-/tt')Aff». 

d .,^ „ ,dH „. / 1 . A dH' 



» /I iT\ TT , i dfi „, / 1 , h dH\ 



= fi' {1 - 2M + A.' + A (/*- A)} 
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= 3 (1 - fih) E' (/i - h) -fiH'. 

= - S/i-ff* + 3 {/* (1 + A*) - 2/t'A} fl* 
= - 3/t {fi» - (1 - 2iih + h') H'} 
= 0, since 1 - 2/iA + A' = 5^. 

Therefore | {d " /*') f } + /^ d^. W = 0. 

This may also be shewn as follows. 

If X, y, z be the co-ordinates of any point, / the distance 
of a fixed point, situated on the axis of z, from the origin, 
and It be the distance between these points, we know that. 



and that 



■■©=«• 



Now, transform these expressions to polar co-ordinates, 
by writing 

a: = rsin^cos0, y = rsin^sin^, ;g:=:rcos^, 

and we get 

and the differential equation becomes 
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or, putting cos ^ = ft. 

Now, putting r = z'h, we see that 

S^ 1 

1 H 

5 = - 



cP /r\ h €P 



.'. the above equation becomes 

4. Having established this proposition, we may proceed 
follows: 

K Pi be the coefficient of A* in the expansion of H^ 

.'. hH=h-^pJb^+pJi!' + ...+Pthf^^ + ..^ 
.\ h^,{hH)^l.2p,h + 2.SpJi^+...+i(i + l)pJi'+... 

Also, the coefficient of h* in the expansion of 
Hence equating to zero the coefficient of h\ 
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Also PtissL nkdoDal integral fanclioa of /k 

And, when M = l, JT= U - 2i -h A V^ 

= 1 +i-h **+.-.+*' + .. . 

Or wten/* = l,/>4 = l. 
Therefore j>^ is what we h&Te already denoted by P^ 

We have thus shewn that, if 4 he leas than 1, 

If & be greater than 1, this soies becomes diTergent. 
Bat we may write 






since ^ is less than 1, 



P, P, P, 

Hence P^ is also the coefficient of i'~^*'^^ in the expan- 
sion <^ {1 — ^^-f i^ ^ in ascending powers of ^ when 4 is 

greater than 1. We may express this in a notation which is 

strictly continuous^ by ;>ayiDg that 

This might hare been anticipated, firom the fact that the 
fundamental differential equation for P4 is unaltered if 
— ^1 + 1) be written in jdace of 1; for the only way in 
which I aj^ars in that equation is in the coefficient of Pp 
which is i \i 4- 1). Writing — [i + 1) in place of i^ this he- 
comes — (i + 1; [— (» + !} + 1; or (i+l) i, and is therefore 
unaltered. 
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5. We shall next prove that 

P-( 1V_^ ''' c^\ 

. '"^ ^' 1.2...td?W' 
where f' = a? + i^+^. 

Let -J. = (*• + y* + «•)■* =/(^), 

and let k be any quantity less than n 

Then {a? + y^ + (z-^ jk)'}"* =f{z^k), 

and, developing by Taylor's Theorem, the coeflScient of k' is 

(-•)'i^r"(-»'rir*(;)- 

Also {a? + y'+iz-ky}-^ = {T'-2kz + k*)-<' 

^n the expansion of which, the coefficient of k* is 

P. 

Equating these results, we get 

-^^ e? /l 






The value of P, might be calculated, either by expanding 
(1 - 2fih + A*)"^ by the Binomial Theorem, or by effecting the 

differentiations in the expression (— 1)' ^ ^ •:7"i(-) > 

L , Z , o ,,, taz \r/ 

mm 

^d in the result putting - = /a. Both these methods how- 

?^cr would be somewhat laborious ; we proceed therefore to 
^^estigate more convenijent expressions. 
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6. The first process shews, by the aid of Lagrange*: 
Theorem, that 

1 //* 

Let y denote a quantity, such that 

h, being less than 1> 

Then 

1 

dy h 1 

Hence, by Lagrange's Theorem, 

therefore, differentiating ■with respect to ft and observing 
that 

^=(l-2M + A')-i, 

(i-2M.A«)-^=i..|(^).^,|.(^-7.... 

+ i.2...iV/t^ 2 ;+•"• 
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7. From this form of P, it may be readily shewn that 
tlie values of yx, which satisfy the equation P^ = 0, are all real, 
ajad all lie between — 1 and 1. 

For the equation 

(/x*— iy = has t roots =1, and t roots = — 1, 

•'- -T- (/x*— 1)*= has i — 1 roots = 1, (i — 1) roots = — 1, and 
one root = 0, 



^— i(/ir*— 1)' = has (t — 2) roots = 1, one root between 1 

ajad 0, one between and = — 1, and (i — 2) roots = —1, 
ctxid 80 on. Hence it follows that 

^~~i {ji* — 1)*= has ^ roots between 1 and 0, and - roots be- 
tween and — 1, if t be even, 

* 1 * 1 

^uid —5— roots between 1 and 0, — ^— roots between and 

1, and one root = 0, if i be odd. 

It is hardly necessary to observe that the positive roots of 
Gach of these equations are severally equal in absolute mag- 
^tude to the negative roots. 

8. We may take this opportunity of introducing an im- 
j Portant theorem, due to Rodrigues, properly belonging to 
I the Differential Calculus, but which is of great use in this 
[ subject. 

The theorem in question is as follows: 
Ifm he any integer less than i, 

^<''-')'° l.i:f;:j »'-'rg.(^-^i-- 
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It may be proved in the following manner. 

If (a?*— 1)* be differentiated t -m times, then, since the 
equation 

(x»-l)' = 

has { roots each equal to 1, and % roots each equal = — 1, it 
follows that the equation 

%«(^-iy=o 



dx 

has i — (i — m) roots (i. e. m) roots each = 1, and m roots 
each = — 1, in other words that {a? — 1 ;"* is a factor of 

We proceed to calculate the other factor. 

For this purpose consider the expression 

(a? + aj (x-^a^ ... (a; + a,) (x + ^SJ (x + fi;) ... (x + fi,). 

Conceive this differentiated (I) i — m times, (II) i + m 
times. The two expressiona thus obtained will consist of an 
equal number of terms, and to any term in (I) will corre- 
spond one term in (II), such that their product will be 
{x + a^ (x + a^) ... (^+a,) {x + 0^) (x + jS^) ... (x + ^j),i.e. the 
term m (II) is the product of all the factors omitted from 
the corresponding term in (I) and of those factors only. 
Two such terms may be said to be complementary to each 
other. 

Now, conceive a term in (II) the product of p factors of 
the form x + a, say a? + a', x + ol\..x + cl^\ and of q factors 
of the form x + fi, say x + ^^^ x + ^^^...x + fi^^y We must 
have^H- j = i — 7n. 

The complementary term in (I) will involve 

p factors x + ^,x + ^\..x + ^\ 

q factors x + a^, x + a^, ... x -V a^,i\. 
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Now, every term in (I) is of i + m dimensions. We have 
accounted for p+ j (or i — m) factors in the particular term 
we are considerinor. There remain therefore 2m factors to 
be accounted for. None of the letters ) 

a, a' ... a , ^,, ^^,,..fi^g^ 

can appear there. Hence the remaining factor must involve 
wt a's and m /8's, — say, 

jOf, gOt , , , ^If 
lH> 2^ • • • wP« 

There will be another term in (II) containing 

{x + ^) (^ + i9") ... (x + ^O {x + a) (x + a,) ... (x + aj. 

The corresponding term in (I) will be, as shewn above, 

(x + a') (x + a") ... {x + u^) (x + ^) (x + y8J ... {x + ^J 

{x + ^a)(x + ,oi)...{x'{r^'x){x + ^/3){x4-j3)...{x + j3). 

} Hence, the sum of these two terms of (I) divided by the 
sum of the complementary two terms of (II) is 

' (« + !«) (« + ««) - (^ + m«) (a? + i/3) (x + ^fi) ... (a?+ JS). 

•r Now, let each of the a's be equal to 1, and each of the /3's 
equal to — 1, then this becomes (^ — 1)"*. The same factor 
euters into every such pair of the terms of (I). Hence 

^^ dx'-^ ^ =(a?'-l)*" ~^, , to a numerical 

fector pris. 

The factor may easily be calculated, by considering that 
the coefficient of x'"^ in ^^^^^^^is2i(2i-l)...(i+m+l), 

ftnd that the coefficient of x*"^ in ^r-r^ — - is 

2i(2i^l)...(t + m+l){i + m)...{,i-rti\\\. 
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Hence the factor is 



1 1.2... (i — m) 

, or ^ ^ 



(i + m) (i + m — 1) ... (i — m + 1)' 1.2... (t + w) ' 

9. This theorem affords a direct proof that C-^^ (/li*—!)*, 

C being any constant, is a value of / (/i) which satisfies the 
equation 

=t(t+i)^.(^*-i)' 

from above, 

■•• |[''''-«i{sS(^"-'>'}]-'(*'+'HI''^*-'''}- ' 

or 

I [o - "•) I {|. <^' - i)j] +••(•>') {I (^' - » •} - «■ 

Hence, the given differential equation is satisfied by put- 

Introducing the condition that P, is that value of f (ji) 
which is equal to 1, when /x = 1, we get 

1 d* J , 

-^*'2*.i.2...;5^^^'"^^'- 

10. We shall now establish two very important proper^ 
ties of the function P, ; and apply them to obtain the develop- 

J22cut of I) in Q. series. 
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The properties in question are as follows : 
If land Txi he uneqiuil positive integers, 

And \[p^d^ = ^ . :^ --:, 

The following is a proof of the first property. 
We have 

Multiplying the first of these equations by i^ the second 
by P^ subtracting and integrating, we get 

Hence, transforming the first two integrals by integration 
by parts, and remarking that 

i(i+l) — m(m + l) = (i — m) (i+ m + 1), 
we get 

<• -"•> ('-f -^. t) -/(> -"•' (tf -f f )^'' 

+ (i - m) (i + m + 1) jPtPJ/i = 0, 



or 



(l-M»)(p.^'-P.^») + (»-m)(i+m + l)Jp^,d'/*=0. 

since the second term vanishes identicaUy. 

r.s. ^ 
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/: 



Hence, taking the integral between the limits — 1 and 
+ 1, we remark that the factor 1 — /t* vanishes at both limits, 
and therefore, except when i — m, or i+fn+ 1 = 0, 

-1 
We may remark also that wo have in general 

pdP,_p dP„ 

a result which will be useful hereafter, 

11, We will now consider the cases in which 

i — m, or i + m + 1 = 0. 

We see that i + m 4- 1 cannot be equal to 0, if » and m are 
both positive integers. Hence we need only discuss the 
case in which m = i. We may remark, however, that since 

Pi == P_(,+i), the determination of the value of 1 PJ cZ/a will also 

fi ^" 

give the value of I F,P_^^^^^ dfi. 

The value of I P*dfi may be calculated as follows : 

(l-2fih+h*)-^ = P,+P,h + ... + PJi' + ...; 

.-. {l-2fih + hY = {P, + PM-' +Pih* + ...y 

= Po' + P^'h* + ... + P,Vi"+ ... 
+ 2P,PJi + 2P,PJi* +... + 2P^PJi' + ... 

Integrate both sides with respect to /a ; then since 



/ 



{l-2fih + hYdfi>=-^logCL-2fih + h'^, 



2h 
we get, taking this integral between the limits — 1 and + 1, 

all the other terms vanishing, Yiy l\ie \jQ60Tcai ya&t -ijtoved. 
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Now log^ = 2(A+|*+... + ^ + ...). 
Hence 2(1 +^ + ... + 2^ + ...) 

Henc^ equating coefficients of A", 



fi 2 

12. From the equation I PjP„j(i/i = 0, combined with 

the fact that, when /t = 1, P| = l^ and that P, is a rational 
integral function of fi, of the degree i, P, may be expressed 
in a series by the following method. 

"We may observe in the first place that, if m be any 

integer less than i, I ii^P^dfi = 0. 

For as P«, P^.^ . . . may all be expressed as rational in- 
tegral functions of ft, of the degrees m, w — 1 ... respectively, 
it follows that yT will be a linear function of P^ and zonal 
harmonics of lower orders, /a"*"* of P^j and zonal harmonics of 

lower orders, and so on. Hence I/i*"P,(//a will be the sum of 
a series of multiples of quantities of the form I FJPd^^ 
m being less than j, and therefore I /x.'"Pj(7/i = 0, if m be any 

J ""X 

integer less than i. 

Again, since 

{l-2nh + h')-i = P, + PJi+ ... + P,h' + ... 

it follows, writing — h for h, that 
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And writing — /a for /a in the first equation, 

P/, Pj'...P/,,. denoting the values which P^, P^...Pi, 
respectively assume, when — /i is written for /a. Hence 
P^^^Pi or — P,, according as i is even or odd. That is, 
P, involves only odd, or only even, powers of ^r,. according r 
as i is odd or even*. 

Assume then 

P, = il,/i* + ^^/A*^+... 

Our object is to determine -4,, -4,,, .... 

Then, multiplying successively by /i'"", /i*^, ... and inte- 
grating from — 1 to + 1, we obtain the following system of 
equations : 



T ~^» 5" T • • • T 'cT^ o T T • • • ~~ "> 



2i-l • 2i-3 2i-25-l 

-4< , -4,_, -4<_5, 



I >v • ^ T • . . "t* gT"; g) Q "t* • • • "^ ^> 



2t_3 ' 2i-5 2i-2s-3 



■"« ■ -^j-j J j_ -'^j-ii 



T" ^ • c% o *• • • • "•" C\ ' A ■% *• • • • ~~ ^* 



2i-2s-l ' 2i-25-3^ ' 2i-45-l 



And lastly, since P^ = 1, when yx = 1, 

^, + ^iHi + . .. + ^,.3. + ... = 1 ; 
the last terms of the first members of these several equa- 
tions being 

AAA 
j-^, -r-^...--^,^„ifibeeven, 

AAA 

JZ2' 5:^4 — y'^A'^^^^^^^- 
is. The mode of solving the class of systems of equa- 
tions to which this system belongs will be best seen by 
considering a particular example. 

* This is also eyident, from the fact that Pi is a constant multiple of 
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Suppose then that we have 



CO V z ^ 



a+a 6+ a c+a 



a + iS 6 + )8 c + ^ 

a + o) ft + o) c + fi> (» ' 

From this system of equations we deduce the following, 
6 being any quantity whatever, 

X y z ^ 1 (g-tt)(g~/9)(a+Q>)(&fQ>)(c+o)) 

a+tf J + dc + ^ oi (w-a)((»-^)(a+^(6+^)(c+5)' 

For this expression is of — 1 dimension in a, 6, e, a, /S, 7, 
^, o); it vanishes when ^ = a, or ^ = )8, and for no other 

finite value of 6, and it becomes = — , when ^ = o>. 
We hence obtain 

rlfrr^^^f y f ^ U ^ (^-a) (g-^ (a+ft)) (6 + <o) (c+O) ) 

and therefor-e, putting ^ = — a, 

1 (a + a)(a + )8) (a + o)) (6 + ft») (c + o) 
G>(a — 6)(a — c) (® — a)(ft)— /8) * 

with similar values for y and 2;. 

And, if <» be infinitely great, in which case the last 
equation assumes the form x +y + z=l,we have 

^ (a + a) (a + 13) 

with similar values for y and «• ^ 

14. Now consider the general system 

• ^— +-%-+...+-%- + ...=0, 



0/+<t4 <%L« + *f4 .««-»•+«»-» 
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""* - ^'-. +...+ ^'-^ +...^0, 



<*< + ««-»• ««-« + a<-it «i-a.+ a<-«. 



^^ +-^^=^+...+ ^*-r +...=- 



a^ + o) a^^ + m (^i^ + ^ ® 

the number of equations, and therefore of letters of the 

forms X and a, being —5— if i be odd, ^ + 1 if t be even ; and 

I — 1 . 

the number of letters of the form a being — ^ if $ be odd, 

and ^ — 1 if i be even. 
We obtain, as before, 



^«-8»~" 



_ j^ (g-g0(g-O ...(g-g,,J ... (a,+tt))(a^4-a))...(g^4-fi))... 
ck> (o)-a5"(a)-aJ)...(a)-a^J... (a, +^)(«i-2+^) •••(««->. + ^)--- 
lind, multiplying by a^_^ + ^, and then putting ^ = — a^^, 

0) (a)-ai) (fi)-a^^...(o>-a,^)... 

15. To apply this to the case of zonal harmonics, we see, 
by comparing the equations for x with the equations for A^ 
that we must suppose o» = x ; and 

a, = t — 1, a^ = i — 3,...a4.^ = i — 2«— !.•• 
Hence 

'-«•'" (-a?) {-(a?-2)}...{(i-2«-l) or (i-2«)} 

^. ^., (2i-2^-l)(2i-2g~3)...{2(t-28)-l]... 
^■"^&(2^-2).,.^x2;,4i..\^-^-\^^x(^~2«)• 
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Or, generally, if f be odd, 

. (2t-l)( 2t-3)...(t + 2) 
' 2.4...(i-l) • » 

J. (2t~3)(2t-5)...t 

*-* 2.4...(t-3)x2 ' 

. _ (2t-5)(2f-7)...(t-2) 
*-* 2.4...(t-5)x2.4 ' 



^ _ , ^ .-J- t(t — 2)...3 
^ "^ ^ 2.4...(i-l)- 

And, if i be even, 

^ _ (2e-l)(2f-3)...(t+l) 
' 2.4...i 

A (2t-3)(2t-5)...(t'-l) 

'-• 2.4...(i-2)x2 ' 

. _ (2^•-5)(2^•-7)...(^•-3) 
'-« 2.4...(i-4)x2.4 ' 



. _ i (»--l)(t-3)...l 

^•~^~-^^ 2~i:::i 

We give the values of the several zonal harmonics, from 
P, to P„ inclusive, calculated by this formula, 

P. = l. 

-^» 2 '^ 2 ' 

_V-i 

~ 2 ' 
p_5 ,_§ 

* 2^ 2'* 
5^» - 3/t 
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^'"2.4.'* 2:2'*^2:i, 

Sofi*-SOfi* + S 
~ 8 ' 

' 2.4^ 2.5J^ ■^2.4'^ 

_ 63/t'-70At'+15/t 

8 ' ' 

„ 11.9.7 , 9.7.5 4 7^5.3 . 5 .3.1 
'"2.4.6'* 2.4x:i'* ■*'2x2.4'* 2.4.« 

_ 23 V- 315/t*+ lOofi*- 5 

16 

p 13.11.9 , 11.9.7 5..9^ij5 , 7.5.3 
' 2.4.6 '^ 2.4 x2'* ■'"2x2.4'* "2.4.6'* 

_ 429/ - 693/*' + 315ft» - 35/t 

16 

15.1 3.11.9 g 13.11.9.7 ^ 11.9.7.5 .- 
' 2.4.6.8 '*"* 2.4.6x2'* ■^2.4x2.4'* 

9.7.5.3 , ^ 7.5.3.1 
2x2.4.6'* ' 2.4.6.8 

6435/1" - 12012/ + 6 930/t* - 1260/** + 35 

f28 * 

„ 17.15.1 3.11 , 15.13.11.9 T . 13jllj.9^ s 
'~ 2.4.6.8 '* , 2.4.6x2 '* ■'■2.4x2.4 '* 

11.9.7.5 , 9.7.5.3 
2x2.4.6'* "'■2.4.6.8'* 

_ 12155/- 25740/ + 18018/- 4620/ + 315/* 

128 ' 

p ^1 9.17.15.13 .11 „ 17.15.13.11 .9 , 1 5.13.11.9.7 , 
" £.4.6.8.10 '* "2.4!.6.%xt'*'^^A.^x^.4i'* 
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13.11:9.7.5 - 11.9.7.5.3 , 9.7.5.3.1. 



/^"^Ov^o A. a q/^ 



2.4x2.4.6'^ '2x2.4.6.8'^ 2.4.6.8.10 

_ 46189/i^"- 109395/i^+ 90090fi'- 500S0fi'+ 3465/^'- 63 

256 

It will be observed that, when these fractions are reduced 
to their lowest terms, the denominators are in all cases 
powers of 2, the other factors being cancelled by correspond- 
ing factors in the numerator. The power of 2, in the 
denominator of P^, is that which enters as a factor into the 
continued product 1 . 2...i. 

16. We have seen that / fi'^F^.dfju — O, it m be any 
integer less than t. 

It will easily be seen that if m + i be an odd numb w, the 
values of j/i^P^.dfju are the same, whether [m be put = l or 

— 1 ; but if w + » be an even number, the values of I/a*" P^ . dfjL 

corresponding' to these limits are equal and opposite. Hence, 
(m + i being even) 

P fjrP,.dfi=^2pfi'^P,.dfjL, 

and then I fjJ^P^.dfi^Of if m = i — 2, i — 4 

Jo 



We may proceed to investigate the value of / /i*" P, . d/i, 

if m have any other value. For this purpose, resuming the 
notation of the equations of Art. 13, we see that, putting 

s m + 1, and © = x , we have 

• - ...... 

5« -1- ?^-? + + ?*=?! + 



ai + w4-l a,^4-m + l at^ + m + 1 

_ (m + l — a<)(m4-l --a<^) ... (m + 1 — g^.,,) ... ^ 



/, 
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and therefore, putting Xf=Af..., a,^i..., a,ssi—l...f 
we get 

X - A A A 

_ (m-i+2)(m-i+4!)...(m-l) . ., ,, 
~ (». + » + l)(«i + »-l) ... («i + 4)(m + 2) " * "« o°a, 

and (CT-» + 2)(iB-t+4)..-m jf/ be even 

In the particular case in which m = t^ we get 

''''^ " (2i+l)(2f-l).'y(i + 3)(i+l) (* ^^^^)- 

17. We may apply these formulaB to develope any positive 
integral power of /i in a series of zonal harmonics, as we 
proceed to shew. 

Suppose that m is a positive integer, and that fiT is de- 
veloped in such a series, the coefficient of F^ being C7|, so 
that 

then, multiplying both sides of this equation by P, and inte- 
grating between the limits — 1 and 1, all the terms on the 

right-hand side will disappear except I C^F^^dfi^ whiclL will 

2 ' 

become equal to ^^ — ^ G^. 



Hence Ci = ?^ J V"^«<?/^ 



which is equal to 0, if m + 1 be odd. Hence no terms appear 
imless m + i be even. In this case we have 



*=(2i + l^( iTP^oEl*, 
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Hence the formula just investigated gives 
r-rgi I 1^ (w-i4-2)(m-i+4)...(m--l) 
^* ^ "*'^(m+4 + l)(m + i--l)...(m + 4)(m + 2) 

if % be odd^ and 

^ - /0--4. 1% (m-t4-2)(m-t + 4),..7yt 

^"■^ "*■ ^(m + »+l)(m+j-l)...(m + 3)(m+l) 

if » be even. 

Therefore if m be odd, 

m^fa , -|\ 2.4.6 ... (m— 1) p 

A* - [Zm + l) ^2^^ j^ ^2m- 1) ... (m + 4j (wn- 2)^'* "*"••* 

^' (m + 4)(m + 2j »"^m + 2 »* 
If m be even, 

/* -^^»» + A;(2m+l) (2m-l) ... (m + 3) (m + l)"^" "•"••• 

^ (m + 3)(m + l) » m + 1 •• 
Hence, putting for m successively 0, 1, 2 ... 10, we get 

tt' = 5— P+-P=^P+ip 
'^ 5.3 "^3 • 3 » 3 •' 

m'-Tc^P + 'P 
'^ 7 . 5 » ^ 6 ' 

= -P+-P 
5 » 6 *' 

ii««9 ^-^ p +5_i_p +1 p 

^—F+-P +-P 
55 ♦ 7 » 5 •• 



"58 .ZONAL HARMONICS. 

"63 » 9 »'**7 " 

'*' = ^^ 13 . 11 .9 . 7 ^•■*' ^ ir~9:7 ^« ■*■ ^ 977^* ■•" 7 -^^ 
• =i6p+24p 10 1 

' »7.^ig 2.4.6 p , 11 ^'^ fi i7 ^ P I ^P 
^ "^^15.13.11.9^'^^^ 13. 11.9^»^^ 11.9^-^9^'' 

"^ 429 '^39 «^ 33 » 3 *' 

8- 17 2.4.6.8 p TO 4.6.8 p 
'^ -.^Vl7.15.13.11.9^«"^^'^15.13.11.9^* 

. Q _6^_8_ „ -_8_ „ 1 p 
^ 13.11.9 *^ 11.9 »^9" 

_128p _64 .48 .40 1 

~ 6435 * T 495 •^« ^ 143 * ^ 99 » ^ 9 •' 

/.'=19,_l^.^-,! „P, + 15 ^-^-^ " 



19.17.15.13.11 •' 17.15.13.11 '' 

+ 11 6.8 p 8 p . 3J.P 

^ 15.13.11 »^' 13.11 "^11 " 
_ 128 192 p ,16p , 56 3 

12155 "2431^ 6'5»^ 143 "^11^' 

„w _ 21 2.4.6.8 .10 p - 4.6.8. 10 p 

'^ "21.19.17.15.13.11 '•'^19.17.15.13.11 » 

■ ■■o 6.8.10 p 8. 10 p 10 p , 1 p 

"•■^*'l7.15.13.11 '"^^ 15. 13.11 ■^*'*' 13.11 »"^ 11^" 

256 p .i28p 32^p 48p ,iO Ip 



^^159 "^2717 • 187 .* \^- * 14a » ■ 11 
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18. Any zonal harmonic P^ may be expressed in a finite 
series of cosines of multiples of 0, these multiples being 
t^, (1-2)^.... Thus 

(l-2M + A»)-i=P, + P,A + ...+P;i'+...; 

therefore, writing cos for /a, and observing that 

1 - 2 cos ^A + A« = (1 - he^^ *) (1 - Ae-^^^^), 
we obtain 

or 



(l + |/i€>^-l^ + |^A«€V-12^ + ... 



^ 2.4...2i ^^ ^-J 

2.4... 2i / 

whence, equating coefficients of U, 
„ 1.3...(2i-l)„ .. , 1.3...(2t-3)l„ ,. _,. 

/•= •^:47:2r^2cosz^+2:^-|^-2cos(t-2)d 

the last term being \ - a 1 — = — r if » be even, and 
1^3^^^(i + l) 1 . 3 ... (i-2) „ a 'fx. j i 

19. Let us next proceed to investigate the value of 

cos mO sin d0. 



f"-' 

J9 
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This might be done, by direct integration, from the above 
expression. Or we may proceed as follows. 

The above value of P, when multiplied by cos m6 sin 6 

(that is by 5 {sin (m+ 1) ^ — sin (m — 1) 0]) will consist of a 

series of sines of angles of the form {% — 2n + (m + 1)} d, that 
is of even or odd multiples of ^, as i + m is odd or even. 
Therefore, when integrated between the limits and ir it 
will vanish, if t + m be odd. We may therefore limit our^ 
selves to the case in which t + m is even. 

Again, since cos mO can be expressed in a series of powers 
of cos 6, and the highest power involved in such an expression 
is cos 6, it follows that the highest zonal harmonic in the 

development of cos md will be P^. Hence / P, cos mO sin dO 

will be = 0, if m be less than u 

Now, writing 

P,= C;cosi^+ (7^ cos {{-2)0 + ... 

we see that Pf cos mO sin ^0 will consist of a series of sines 
of angles of the forms (m + 1 + 1) ^, (m + i — l) ... down to 
(m — t — 1)^, there being no term involving m0, since the 
coefficient of such a term must be zero. Hence 



Jo 



P, cos m0 sin d0, 



will consist of a series of fractions whose denominators in- 
volve the factors m + i+ 1, w + 1 — 1 ... m—t — 1 respectively. 
Therefore when reduced to a common denominator^ the result 
will involve in its denominator the factor 

(m + t + l)(m+t-l)...(m + l)(m-l)... (m-t-1) 
if 7» be even, and 

(m + t + l)(m + t-.l).-(m + 2)(m-.2).,.(m-t-l) 
if m be odd. 

For the numerator we may observe that since 

P^ coa m 6 sin flcE^ 



/■ 

Jo 
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vanishes if m be less than i, it must involve the factors 
m — (i — 2), m— (i— 4) ... »i + (i — 2), and that it does not 
change sign with nu Hence it will involve the factor 

{m- (t-2)} {m- (i-.4)} ... (m -2) m* (m + 2) ... (m + ?- 2) 
if m be even, and 

if m be odd. 

To determine the factor independent of m, we may pro- 
ceed as follows : 

P, = C;cos {0 + G^ cos (t- 2) ^ + ... ; 
/. JFJcos w^ = 5 C, {cos (m + i) 5 + cos (m — 1)6] 

+ 5 CL, {cos (m + i- 2) ^ + cos (m- i + 2) 5} + ... ; 

.-. P|C0sm^sin^ = T (7,{sin(m + i + l) ^ — sin (m + i-1)^ 

+ sin (m — t +1) ^ — sin (m — i — 1) 5} 

+ sin(w — 4 + 3)^ — sin (m — i+1) 5} + .., ; 
COS mO sin ^ dQ 



•■■ //• 



_c.{^ 1— + -1 l_l 

2 (m + i + 1 m + i — 1 m — i+1 m — t — IJ 

. <5=,|^^ ^— +^^ ^— 1+ . 

^ 2 (m + t-l m + i-3^«t-t>3 w-i + lj 

~ ^' I w" - (i + 1/ ^ TO* - (» - l)'j 
_ f »-l . t-3 1 . 
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Now, when m is very large as compared with t, this be- 
comes 



= -2 



c7< + c7^2 + .., 2 



31 S > 

m m 



since C,+ (7^, + ... = 1, as may be seen by putting ^=0. 

• . r» . 2 

Hence I P, cos m^ sin dO tends to the limit — « , as m 
Jo ^?* 

is indefinitely increased. 

The value of the factor involving m has been shewn 
above to be 

{m-(i - 2)} {m- (^-4)} ... (m-2) m'(m + 2) ... (m +t- 2) 
{m-(^ + l)j{m- (i-l)} ... (w- 1) (m + 1) ... (m + t + l) 
if m be even, and 

{m-(i-2)}{7n-(t-4)}...(m-l)(m + l)... (m + t-2) 
'jm - (i+1)] [m - (i- 1)} ... (w -- 2) (7/1 + 2) ... (m + i+ 1) 

if m be odd. 

Each of these factors contains in its numerator two factors 
less than in its denominator. It approaches, therefore, when 

m is indefinitely increased, to the value —5 « Hence 



Jo 



m 
cos m6 sin d0 





{yn-(i-2)|{m-(/-4)]...(m-2)m'(yn+2)...{m+(i-2)} 
[m-ii-^-Vj] [m-(i-l)j ... (m -- 1} (w + 1) ... [m + (i + 1)} 

if m and i be even, and 

.^ {rn^f-2[} { m-(i-4)}...(m~l)(m+l)...{m+(i-2)} 
{7w-(^•+l)Hm-(V--l)}...(m--2)(m+2)...{m+(^+l)} I 

if m and i be odd. 

In each of these expressions i may be any integer such 
that w ■- 1 is even, i being not greater than m. Hence they 
will always be negative, except vsKeu i xs equoZ to m. 
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20. We may apply these expressions to develop cos mO 
in a series of zonal harmonica. • 

Assume 

coswe = 5«P^+5,^^ + ...+5,P, + .- 

Multiply by P« sin 0, and integrate between the limits 
and TT, and we get 

^ {m - {i- 2)} {m - (^ - 4)} ... {m + {{- 2)} 2_ 

{m - (t + 1)} [m - (i-l)] ... [m + {i+1)} ~ 2t + 1 '' 

Hence 

« rgj , ,>. W-?t-2)}[m-(z-4)}...(m + (t-2)} 

' '^ "^ ''{m-(t + l)j{m-(i-l)} ... {»» + (t + l)} 

Hence, putting m successively = 0, 1, 2, ... 10, 
cosO^ = P,; 
coa$ = P^; 

co8 2^ — 6-^^^P.-|p. 
3 • 3 ^' 

5 » 5 " 



-lJ-7> 
3.5 • 



35 * ^1 » 15 •' 
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^. ^^ 2.4.6.8 ^ ^ 4.6 
cos5g=:~ll . . ^ ^^r-^— T-. -Pk-7 



-1.1.3.7.9.11^ 1.3.7.9 » 



-^sV' 



COS 



65 = -13 



-1?? p_§ p _lp . 
63 » 9 » 7 " 

2. 4. 6'. 8. 10 p 

-1.1.3.5.7.9.11.13 « 



4.6'.8 6' p 1 

1.3.5.7.9.11 V 3.5.7.9 « 5.7 • 

512 384 £ J, 
231 • 385 * 21 « 35 •' 

„a tr 2.4.6.8.10.12 _ 

cos75=-15_^^ 3^ g ^j j3 ^gP, 

_ 4.6.8.10 y 6.8 p 3^ p 

1.3.5.9.11.13 » 3.5.9.11 » 5.9 ' 

1024 128 112 J^ 

■'■1 llT-'s A.CtK-'-t 1K-^I> 



429 ' 117 » 495 8 15 

2.4.6.8M0.12.14 
-1.1.3.5.7.9.11.13.15.17 



cos8g = -17 , -f'T'," •/"•/,•/,, ,^ P8 



4.6.8M0.12 p_ 6 . 8M0 „ 

1.3.5.7.9.11.13.15 ' 3.5.7.9.11.13 • 

K 8' p l_p 

5.7.9.11 » 7.9 • 

16384 p 4096 p 256 p _ j64 p_J^p 
~ 6435 « 3465 ' 1001 * 693 » 63 • ' 

2.4.6.8.10.12.14.16 



cos 95 = - 19 



- 1 . 1 . 3 . 5 . 7 . 11 . 13 . 15 . 17 . 19 » 



4.6.8.10.12.14 p 6.8.10.12 „ 

1.3.5.7.11.13.15.17^ 3.5.7.11.13.15 • 

7 8-10 p 3 1 P 
b.T.W.lS » 7.11 ' 



COS 105 = - 21 
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32768 3072 p 128 16 „ 3 p 
121§5 • 2431 ' 455 * 143 » 77^' 

2.4.6.8,10M2.14.16.18 



-1.1.3.6,7.9.11.13.15.17.19.21 " 
^ 4.6.8J0M2.14.16 

1,3.5.7.9.11.13.15.17.19 » 



6.8.10M2.14 8.10M2 

3.5.7.9.11.13.15.17 ' 5.7.9.11.13.15 * 

7.9.11.13 » 9.11 * 

_ 131072 _ 32768 p _5U J^ p _ ^00, p 

"46189 " 24453 • 1683 • 1001 ' 9009 ' 

_JLp 

99 » 

21. The present will be a convenient opportunity for 
investigating the development of sin 5 in a series of zonal 

harmonics. Since sin d=(l— ;t*)i it will be seen that the 
series must be infinite, and that no zonal harmonic of an odd 
order can enter. Assume then 

8in5= C.P.+ G,P, + ... + aPi+ ... 

i being any even integer. 

Multipljdng by Pj, and integrating with respect to /i 
between the limits — 1 and + 1, we get 



/ 



_P,sin5<f/t = 2^q:j-Ci; 



•'• Gi = — 5 — I -Pi sm 6 afjt 

^^"^^ f'Ptsm*ddd, 
Jo 



2 .,0 

supposing Pi expressed in terms of the cosines of 6 and its 
multiples 

= ?»±I f'p^ (1 _ cos 20) d0. 
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Hence, putting t = 0, 

1 3 

Putting i = 2, and observing that -^2 == 7 + J ^^ 2^> 



4 ' 4 



^ 5 f'^Cl + 3 cos 20) (1 - cos 25) ,^ 
^^^ilo 4 ^^^ 

leio 1 



l + 2cos25-|(l + cos4^lcZ^ 



- 32 '^^ 



For values of i exceeding 2, we observe, that if we write 
for Pi the expression investigated in Art. 18, the only part 

of the expression I P* (1 — cos 20) d0 which does not vanish 

Jo 

will arise either from the terms in P< which involve cos 20, or 
from those which are independent of 0, We have therefore 

* 4 Jo L2.4...(* + 2) 2.4...(i-2) 

_ 2t + l 1.3...(e-l)1.3...(t-3) 
4 •2.4... t 2.4... (»-2) 

*-l . » + l2cos2^)(l-cos2^rf^ 



/, 



\ » i + 2 ^y 



_ 2t + l 1.3... (t-1) 1.3... (i- 3 ) / i-l _ f+1 N 
~ 4 2.4... i 2A...{i-2y\ i i+2j 
2/+ 11.3... (t-1) 1.3... (i-S) 

~ . "^ 2 2A...i{i+2)2A...{i-2)i' 

Hence sin^ = |P,-g|P,- ... 

_ (2t + l)7r 1.3... (i-1) 1.3... (t-3) p _ 

2 2.4...t(i + 2)2.4... (»-2)» * *" 
/ hslng any even integer. 
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dP 

22. It will be seen that -i—* , being a rational and integral 

function of /**"*, /a*"^..., must be expressible in terms of 
-Pi_i> -^«_e»-» To determine this expression, assume 

then multiplying by P«, and integrating with respect to fi 
from — 1 to + 1, 

Now, since i > w. 



r-. 



^^-d'/* = 0; 



dfi 

.-. £ p„ ^' cZ/* = iP^p;\^ - iP„pr = 2, 

, since either w or i must be odd, and therefore either P^ or 
P, =5 — 1, when /L6 = — 1 ; 

2 
.*. 2 = 2^^ (7« or (7,= 2m + l; 

.-. ^' = (2i - 1) P^, + (2i - 5) P.., + (2i - 9) P^, + . . . 

I Hence f-^ = (2-l)i^. 

'; 23. From this equation we deduce 

J It, 

the limits fi and 1 being taken, in order that P, — P^.j may 
be equal to at the superior limit. 
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Now, recurring to the fundamental equation for a zonal 
harmonic, we see that 

• P—P = — ^^~^ (\ — uF\ '^ 



=-^.4T)a-.-)^-&' 



dyi» 



24. We have abeady seen that / T^P^ cZ/a = 0, i and w 
being different positive integers. Suppose now that it is 
required to find the value of I PJP^ dfi. 

We have already seen (Art. 10) that 



And, from aboye, 



^^~'*^"^ — 2m +1 ^^"« '^"-'^* 
CpP da- ^ ( "^(^+l) p(p Ip ^ 
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25. We "will next proceed to give two modes of ex- 
pressing Zonal Harmonics, by means of Definite Integrals. 
The two expressions are as follows : 






{/^±(/A'-l)*coS^r^' 

These we proceed to establish. 
Consider the equation 



1 f » 



d^ 



h cos ^ (a« - J«)4 * 

The only Umitation upon the quantities denoted by a 
and b in this equation is that 6* should not be greater than 
a*. For, if 6* be not greater than a^, cos ^ cannot become 

equal to t while Sr increases from to tt, and therefore the 

expression under the integral sign cannot become infinite. 

Supposing then that we write z for a, and V— 1 p for 6, 
we get 

1 f "• cl^ _ 1 

TtJo «— V~1/)C08^ («*+/>*)*' 

We may remark, in passing, that 
Jo « — V— Ipcos^ h ^ + V — 1/)C0S^ 



Jo 



s«^ » 



z* + p^ cos' ^ 
and is therefore wholly real. 

Supposing thiit /)' = aj" + ^, and that Ci?+y^ + z*^i^, we 
thus obtain 

1 f' d^ ^1 

ttJo 2f-V^pcos^ r' 
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Differentiate i times with respect to z, and there resulU 

r^ (-1)' jl^ f' (?^ 

' ~ TT 1 . 2 . 3...t d«* Jo^ — V— 1 p cos Sr 

TT Jo (z-ATTipCOS^)**** 

In this, write /tr for z, and (1 — /i')* r for p, and we get 

p.ir ^ 

'"7rJo{M-(^*-l)4cos^r^' 
Avhich, writing tt — Sr for ^, gives 

1 r* dI5- 

26. Again, we have 

1 1 r'__d±__ 

(a* — J')* TT Jo a — 6 cos '^' 

In this write l^fih for a, and + (//'* — 1)^ A for b, which is 
admissible for all values of h from up to ^— (/a*— 1)*, and 
we obtain, since a" — 6* becomes 1 — 2fjJi H- h*, 



1 1 r» dyjr 



Acos-^ 



^If^ 

^Jo l-{iE*±(/A'-l)*COS'^} A* 

=i^(^'^[l+{^±(^«-.l)icost}A+... 
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Hence, equating coefficients of h\ 



P,=i J' (^ + (y(i«-l)icosf }*di^. 

The equality of the two expressions thus obtained for P, is 
in harmony with the fact to which attention has already 
been directed, that the value of P^ is unaltered if — (i + 1) 
be written for i. 



27. The equality of the two definite integrals which 
thus present themselves may be illustrated by the following 
geometrical considerations. 




Let be the centre of a circle, radius a, C any point 
within the circle, P CQ any chord drawn through 0, and let 
OC-^h, COP = %COQ=^yfr. Then (7P» = a* + 6* - 2a6 cos S^, 
(7^ = a' + 6^* — 2a6 cos y^. Hence 

(a«+6«-.2a6cos^)(a»+6*-2a5cosi|r) = (a"-6y; 



sinSr^ 



sin -^ flf^ 



••a«4.6«-2a6cos^^a'+6'-2abcoa>|r 



= 0. 
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Again, since the angles OPG, OQG are equal to one 
another, 

sin g_ sin OPC _ sin OQG _ sin '^ 
CF OU " 00 '^ CQ ' 

sin^ sin'^ 

' • (a- +6« - 2ab cos^)* "" (a'* + 6' - 2ab cos ^)* ' 

whence — . i ^ j = 0. 

(a* + 6' - 2a6 cos ^)i ^ (a" + 6' - 2a6 cos f)^ 

^^ / 2 . x2 o ^ ^:^=(a' + i'-2a5costy^-i; 

(a^ + 6*- 2a6 cos S-)*^ ^ 

... (a« - 5y'« , . ,, f\ ^^= - {a' + 5«- 2a6cos t)' df. 

In this, write a' + 6* = /a, 2a& = T (/a* — 1)*, which gives 

£j2 _ J2 _. J^ g^j^^ ^Q g^|. 

da- 
{;.±0.»-l)*cosar = - ^^ ± (^' - ^)* ^^«^J'^^- 

We also see, hy reference to the figure, that as ^ in- 
creases from to TT, -^ diminishes from tt to 0. Hence 

•'O {/i ± (/Lfc* — 1)* cos ar]*+^ JO 

28. From the last definite integral, we may obtain an ex- 
pansion of P, in terms of cos Q and sin 0, Putting /a = cos ^, 
we get 

Pi = 2^;;. J [{cos e + V^ cos -^ sin ^}* 

+ {cos ^ — V— 1 COS -^ sin d}*] d^ 

«- /'/(cos ^)' - '^^J^ coaH C<^^^ 
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TVT f'f ,\fmj, (2m-l) (2»»-3)...l 

And » (t- 1)— (t- 2m + 1) (2m - 1) (2to-3)...1 

1.2..,2»» 2OT(2m-2)...2 

_ i {{-!). ..(i-2m + l) 
(2.4....2m)» ' 

.-. P, = (cos ey - ^^^^ (cos ^" (sin ^« + . . . 

+ (-ir^^^^^§7^^5^^^(cos^)'-(sin^-+... 



CHAPTER III. 

APPLICATION OP ZONAL HARMONICS TO THE THEORY OF 
ATTRACTION. REPRESENTATION OF DISCONTINUOUS 

FUNCTIONS BY SERIES OF ZONAL HARMONICS. 

1. We shall, in this chapter, give some applications of 
Zonal Harmonics to the determination of the potential of a 
solid of revolution, symmetrical about an axis. When the 
value of this potential, at every point of the axis, is known, 
we can obtain, by means of these functions, an expression 
for the potential at any point which can be reached from 
the axis without passing through the attracting mass. 

The simplest case of this kind is that in which the 
attracting mass is an uniform circular wire, of indefinitely 
small transverse section. 

Let c be the radius of such a wire, p its density, k its 
transverse section. Then its mass, M, will be equal to 27rpcfc, 
and if its centre be taken as the origin, its potential at any 

M 

point of its axis, distant z from its centre, will be — r. 

Now, this expression may be developed into either of the 
following series : 

ML lz\1^3z^_ , .. 1.3...(2t-l) g" \ 
Tj 2c* '2.4c* '"^^ ' 2.4...2i c""*""'}*"^^^' 

We must employ the series (1) if a be less than c, or 
if the attracted point lie within the sphere of which the ring 
is a great circle, and the series (2) if « be greater than c, 
or if the attracted point lie mt\io\x\. xJcis s^et^. 
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Now, take any point whose distance from the centre is r, 
ad let the inclination of this distance to the axis of the 
ng be 0, In accordance with the notation already em- 
loyed, let cos 0^ fi. Then, the potential at this point will 
e given by one of the following series : 

I { 1 r^ 1 3 r* 

i JP — _ P 1-4- ■ ' P - — 



,|^o 2^«r^^2.4 V* ••• 

/ -.m I* 3« 5...(2i — 1) p c_ 1 ,^,v 
^^ ^ 2.4.6...2i ^^^^^^--j--^^^ 

For each of these expressions, when substituted for V, 
atisfies the equation V*F=0, and they become respectively 
qual to (1) and (2) when is put = 0, and consequently 
= z. The expression (2') also vanishes when r is infinitely 
reat, and must therefore be employed for values of r greater 
lan c, while (!') becomes equal to (2') when r = c, and will 
lerefore denote the required potential for all values of r 
88 than c. 

These expressions may be reduced to other forms by 
Bans of the expressions investigated in Chap. 2, Art. 25, viz. 

P, = -["(/* + V/^^ cos ^yd^, 

or P,= - r(ji'\^J]F^ cos fy^'^df. 

Substitute the first of these in (1') and (observing that 
= z) we see that it assumes the form 

/•' r l {^+(g*-O^c08^}' 

■*'2.4 c* '"[ ' 
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which is equivalent to 

Mr ^ 

The substitution of the last form of P< in the series [H) 
brings it into the form 



Mfi. 



c^ 



1.3 0^ 

■*'2.4{;3+(«^-r*)icosS^}' 
which is equivalent to 

Mr ^ 



•A 



IT j( 



2. Suppose next that the attracting mass is a hollow shell 
of uniform density, whose exterior and interior bounding 
surfaces are both surfaces of revolution, their common axis 
being the axis of z. Let the origin be taken within the 
interior bounding surface ; and suppose the potential, at any 
point of the axis within this surface, to be 

Then the potential at any point lying within the inner 
bounding surface will be 

AJP^ + A^P^r + A^Py+...+A^P/ + ... 

For this expression, when substituted for F, satisfies the 
equation V^F=0; it also agrees with the given value of 
the potential for every point of the axis, lying within the 
inner bounding surface, and does not become infinite at any 
point within that surface. 

Again, suppose the potential at any point of the axis 
without the outer bounding surface to be 
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Then the potential at any point lying without the outer 
boimding surface will be 

^qP^ I ^1^1 , A-P2 , , -^<-P< . 

r r^ r^ "i" •••■*• ^<+i -r .... 

For this expression, when substituted for V, satisfies the 
equation V F= ; it also agrees with the given value of the 
potential for every point of the axis, lying without the outer 
bounding surface, and it does not become infinite at any 
point within that surface. 

By the introduction of the expressions for zonal har- 
monics in the form of definite integrals, it will be found that 
if the value of either of these potentials for any point in the 
axis be denoted by ^ {z), the corresponding value for any 
other point, which can be reached without passing through 
any portion of the attracting mass, will be 



i r0{« + (2^^ -7^)4 cos ^}d^. 



3. We may next shew how to obtain, in terms of a series of 
zonal harmonics, an expression for the solid angle subtended 
hy a circle at any point. We must first prove the following 
theorem. 

The solid angle, svhtended hy a closed plane curve at any 
point, is proportional to the component attraction perpendicular 
fe the plane of the curve, exercised upon the point by a lamina, 
of uniform density and thickness, bounded by the closed plane 
curve. 

For, if dS be any element of such a lamina, r its distance 
from the attracted point, the inclination of r to the line 
)€rpendicular to the plane of the lamina, the elementary 
olid angle subtended by d8 at the point will be 

dS cos 

And the component attraction of the element of the 
imina corresponding to dS in the direction perpendicular 
its plane will be 

yoifc-^COS^, 
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p being the density of the lamina, Jc its thickness. Hence, 
for this element, the component attraction is to the solid 
angle as pk to 1, and the same relation holding for every 
element of the lamina, we see that the component attraction 
of the whole lamina is to the solid angle subtended by the 
whole curve as pk to 1, 

Now, if the plane of ojy be taken parallel to the plane 

of the lamina, and V be the potential of the lamina^ its 

component attraction perpendicular to its plane will be 

dV 
— -i—. Now since Fis a potential we have V'F=0, whence 

_V*F=0, or VM-^j=0. Hence --j- is itself a potential, 

and satisfies all the analytical conditions to which a potential 
is subject. It follows that, if the solid angle subtended by 
a closed plane curve at any point (x, t/, z) be denoted by 
o), G) will be a function of x, y, z, satisfying the equation 
V^o) = 0. Hence, if the closed plane curve be a circle it 
follows that the magnitude of the solid angle which it sub- 
tends at any point may be obtained by first determining 
the solid angle which it subtends at any point of a line 
drawn through its centre perpendicular to its plane, and 
then deducing the general expression by the employment 
of zonal harmonics. 

Now let be the centre of the circle, Q any point on the 
line drawn through perpendicular to the plane of the 
circle, JS any point in the circumference of the circle. With 
centre Q, and radius QO, describe a circle, cutting QE in L 
From L draw LN, perpendicular to Q 0. 

Let OE^c, OQ = z. 

Then^i = (c^ + ^)i.^, OA'^^^rf^ {(c' + ^)i-;.} 



{^\%*]i' 
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And the solid angle subtended by the circle at 

ON 



= 4inr 



2^ 



= 27r(l ^—... 



To obtain the general expression for the solid angle sub- 
inded at any point, distant r from the centre, we first 
evelope this expression in a converging series, proceeding 
y powers of z. This will be 

c:. 2c» 2.4c» ^ ' 2.4...2i c"*'"*''" 



:ir|l- 

f ^ be less than c, and 

^^\2z^ ^Az'^'" ^ -"^ 2.4...2i ^'"*"-- 
f 2 be greater than c. 

Hence, by similar reasoning to that already employed, 
ve get, for the solid angle subtended at a point distant r 
rom the centre, 

•'"r" c +2 c' 2.4. c' +••• 

_ /_ 1 ^* 1- 3 ... (2t — 1) -f;;4.i ?• 
*■ -^ 2.4...2i c"** ■*■ 
f r be less than c, and 

(2 r* 2.4 r* ■^•■* '^ ' 2.4...2i r« "*"•••] 
f r be greater than c. 

4. We may deduce from this, expressions for the potential 
f a circular lamina, of uniform thickness and density, at an 
ictemal point. For we see that, if V be the potential of 
ich a lamina, k its thickness, and p its density, we have for 
point on the axis, 



d3 

r.H. 






(c* + z^)^ 



m ' 



50 APPLICATION OF ZONAL HAEMONICS 

whence F= 27rpk {(c* + s^)^ — z} 

if M be the mass of the lamina. 

Expanding this in a converging series, we get 



. iff 



,6 



^■*'2c 2.4c»'^2.4.6c* 



^ -' 2.4.6...2i c"-'"^ 



if z be less than c, and 

3 c' 
6z' 



?t22 2.4«*'^2.4. 



, ,,, 1.1.3...(2i-.S) ^*' 



2.4. 6.. .2i s^' 
if s be greater than c. 

Hence we obtain the following expressions for the 
tential of an uniform circular lamina at a point distai 
from the centre of the lamina : 



4({p..-p. 



^ipy IIP/, 

^2 2.4 c» ^•" 

., 1.1.3...(2i-3)P^r^ , 
^ ^ 2.4.6...2t c«-» "^ 



if r be less than c, and 

c»l2 r 2.4 r" ''"2.4.6 »^ *•' 






2.4. 6. ..2i r^'- 
if r be greater than c. 
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Y" *ttay be shewn that the solid angle may be expressed 
avtlie fori ^ ^ ^ 

JO [c'+{z+(z*-r')i COS 0}li 
*M the potential of the lamina in the form 



jlf^lf^+iz + iz'-r^j^cosOn^de-^. 



5. As another example, let it be required to determine 
'^e potential of a solid sphere, whose density varies inversely 
s the fifth power, of the distance from a given external point 
^ at any point of its mass. 

It is proved by the method of inversion (see Thomson 

Qd Tait's Natural Philosophy, Vol. 1, Art. 518) that the 

M 
otential at any external point P' will be equal to - , p; , 0' 

eing the image of in the surface of the sphere, and M 
16 mass of the sphere. We shall avail ourselves of this 
isult to determine the potential at a given internal point. 

Let G be the centre of the sphere, the given external 
oint. Join GO, and let it cut the surface of the sphere in A, 
ad in GA take a point 0', such that CO.GO-^ GA\ Then 
^ is the image of 0. 

Let P be any point in the body of the sphere, then we 
ish to find the potential of the sphere at P. 

Take as pole, and OG as prime radius, let OF = r, 
Oa=^. Alsolet(7u4 = a, CO = c. 

Let the density of the sphere at its centre be p, then its 
3nsity at P will be /> -^ . Hence 



J/=27r/L^r'sin0drdfl, 



V-'i 



\ 
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the limits of r being the two values of r which satisfy tlie 
equation of the surface of the sphere, viz. 

r* + c' — 2cr cos ^ = a', 

and those of beinff and sin"* - . 

° c 

Hence, if r^, r, be the two limiting values of r, we Lave 

^j 11 2c cos ^ /I 1\ 
Now -^ ^=-j^ ^ 1. 

., 1 . 1 2c cos ^ 
Also — 4- __ s= — ^ _ ^ 



r^ r^ c — a 



r^r^ <? — 0?^ 



. /I 1\ ^ {c'cos'g-(c'-a')}i 
^r, r,/ c* — a* 



9 2 

c — a 



.'. ilf= - ^ -^j 2 . -| g I cos ^ sm 5 (a*— c sm'^)- U 

^ C— tt C ^^ d J (^ 
a 

= A-^^2 r'"'" cos e sin ^ (a» - c» sin» ^ ^ <?5 



(c--a') JO 

4 Trpc* 
~ 3 (c^ - a") 



a\ 



Now, if F be the potential at P, we have (see Chap 
Art. 1) 



rf^(rF) . 1_ d 
dr' ^8in0cl0 
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This is satisfied by F= — ^ — j-. 

Assume then^ as the complete solution of the equation. 



r=-^ 



2irpc^ 



3 r* 



+ (A+^)i^«+(A-+5)i\ + ... 



+ (J/ + ^|^,)P, + .... 



It remains to determine the coefficients -4^, A^.,.A^.,,B^, 

2?j...-B„ so that this expression may not become infinite for 

; any value of r corresponding to a point within the sphere, 

and that at any point P on the surface of the sphere it may 

be equal to jy-pf where CP: OP :: a: c, and therefore, at 

j the surface, 

Tr__ Mc 1 __ 4 irptfa* 

VaP"3(c*-ayr' 

And, at the surface, we have 
r*- 2crcos^ + c'-a' = 0; 

1 ^ 1 /I 2ccosg \ 

. 4 ^pc'q' _2 ^pc* /P. PA / ^\ 

(7? \ 
A^r+ -^M Pj+ ... identically, 

4 TrpcV _/2 TTpc' \p 

3(c'-ay"V3c'-a" V '' 

and ^, ^^...-^...-^^, A^...A^ all = 0. 
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Hence since P, = 1, 






whence we obtain, as the expression for the potential at asy 
internal point, I 

6. We shall next proceed to establish the proposition that 
if the density of a spherical shell, of indefinitely small thid- 
ness, be a zonal surface harmoniCy its potential at any internal 
point mil he proportional to the corresponding solid hx- 
monic of positive degree, and its potential at any external 
point will be proportional to the corresponding solid Jiarmonic 
of negative degree. 

Take the centre of the sphere as origin, and the axis of 
the system of zonal harmonics as the axis of z. Let 6 be the 
radius of the sphere, Bb its thickness, U its volume, so that 
U= 47r6'S6. Let CP^ be the density of the sphere, P, being 
the zonal surface harmonic of the degree i, and G any con- 
stant. 

Draw two planes cutting the sphere perpendicular to the 
axis of z, at distances from the centre equal to f, f+df 
respectively. The volume of the strip of the sphere inter- 

cepted between these planes will be ^^ U, and its mass will be 
CP,§U. 

Now f = hfi, hence cZf = ^^/^> ^^^ this mass becomes 
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56 the potential of this strip at a point on the axis of z, 
at z from the centre, will be 

CU P, , 

2 (i^ + b" - 2bzfi)^ '^' 

h may be expanded into 

^|(p, + P.^+... + P.^;+....)c?/*if^>5. 

'o obtain the potential of the whole* shell, we must inte- 
j these expressions with respect to fi between the limits 
md +1. Hence by the fundamental property of Zonal 
nonics, proved in Chap. ii. Art. 10, we get for the po- 
al of the whole shell 

-^'; — - 77+1 at an internal point, 

Q . — - -j+j at an external point. 

i'rom these expressions for the potential at a point on 
axis we deduce, by the method of Art. 1 of the present 
)ter, the following expressions for the potential at any 
t whatever : 

CU Pr* 

V^ = ^ — ^ TTPi at an internal point, 

,^ CU pp , . 1 • . 

Kj = -^, — =- —^^ at an external pomt. 

i'rom hence we deduce the following expressions for the 
lal component of the attraction on the point. 

!f ormal component of the attraction on an internal point, 
sured towards the centre of the sphere, 

_ dV,_ i f^jrP/-' 
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Normal component of the attraction on an external point, 
measured towards the sphere, 



^. _ i+l. cu^^' 



j+i 



dr 2i + 1 r 

In the immediate neighbourhood of the sphere, where r is 
indefinitely nearly equal to J, these normal component at- 
tractions become respectively 



p 
rtTT * 

2i + l^^ 6*' 



1+1 jyjj -t^i 



and their difference is therefore 

P. 



CU 



V' 



And writing for U its value, 47r&'^S3, this expression be- 
comes 

47rS6 . CP,. 

Or, the density may be obtained by dividing the alge- 
braic sura of the normal component attractions on two points, 
one external and the other internal, indefinitely near the 
sphere, and situated on the same normal, by 47r x thickness 
of the shell. 

7. It follows from this that if the density of a spherical 
shell be expressed by the series 

Cq, G^, G^... Gi... being any constants, its potential (F,)at 
an internal point will be^ 

\ b ^3 b* ^o b" ^'"^2i + l 6'*' ^ ""I 
and its potential ( F,) at an external point will be 

In the last two Articles, by the word *' density " is meant 
''volume density," i.e. the mass of an indefinitely small 
element of the attracting spVieie, ^ynV^^^Vt^ tl^e volume of 
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ame element. The product of the volume density of 
3lement of the shell, into the thickness of the shell in 
leighbourhood of that element, is called "surface den- 
We see from the above that, if the surface density 
[pressed by the series 

otentials at an internal and an external point will seve- 
be 

[b ^"3 b' ^o b' ^•••^2i + l 6*+i ^-7' 



■"{ 






his variation in surface density may be obtained either 
)mbining a variable volume density with an uniform 
ness, as we have supposed, or by combining a variable 
ness with a uniform volume density, or by varying both 
ness and density. 

. We have seen, in Chap, ii., that any positive integral 
r of fi, and therefore of course any rational integral 
ion of /i, may be expressed by a finite series of zonal 
onics. It follows, therefore, that we can determine the 
itial of a spherical shell, whose density is any rational 
ral function of fi. 

uppose, for instance, we have a shell whose density 
3 as the square of the distance from a diametral plane, 
ag this plane as that of a?y, the density may be ex- 

ed by p/M*, or p ts • We have seen (Chap. ii. Art. 20) 

[ence, by the result of the last Article, the potential 
)e 

p ^ (r + K -hr) a** ^^ internal point. 



58 APPLICATION OF ZONAL HARMONICS 

or, since P/' = ^ "" r^ = -^-^ — , we obtain 
p-^ (5 "^ R — P — ) ^^^ *^® potential at an internal point, 
/o o j- + -^ ( — -3 + — g- J Mor that at an external point. 

9. As an example of the case in which the density is re- 
presented by an infinite series of zonal harmonics, suppose we 
wish to investigate the potential of a spherical shell, whose 
density varies as the distance from a diameter. Taking this 
diameter as the axis of z, the density will be represented 

by p sin $, or p (1 — p?)K We have investigated in Chap. n. 
Art. 21, the expansion of sin in an infinite series of zonal 
liarmonics. Employing this expansion, we shall obtain for 
the potential 

irpjl 

2 b 



2^ 



16^6^ ••• 2.4...t(i-f2)-2.4...(i-2)i *6* "j' 

or 

1 n^l^^l^p^^^ 1.3...(t-l) 1.3...(t-3) p V 1 
2^^|2r 16 V ••• 2.4...i(i+2)-2A..(i-2)i V^^ 'i' 

according as the attracted point is internal or external to the 
spherical shell, i being any even integer. All these expres- 
sions may be obtained in terms of surface density, by writing, 
instead of pU, 47rcV. 

10. We may next proceed to shew how the potential of 
a spherical shell of finite thickness, whose density is any solid 
zonal harmonic, may be determined. Suppose, for instance, 
that we have a shell of external radius a, and internal radius 
a', whose density, at the distance c from the centre, is 

yi Pfi^i h being any line of constant length. 

Dividing the sphere into concentric thin spherical shells, 
of thickness c?c, the potenli^A. oi ^ii^ ot^a q£ these shells, of 
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stdius c, at an internal point distant r from the centre will 

>^ obtained by writing c for h, ^ for C7, 4m<?dc for J7, in 

lae first result of Art 6. This gives 

p 4iTrc*dc PfC*r* 4nr P -p t j 

To obtain the potential of the whole shell, we must inte- 
frate this expression, with respect to c, between the limits 
% and a. Tins gives 

Again, the potential of the shell of radius c, at an external 
•oint, will be 



i + lA' ' 



,8if2 



pW^fcP.c-^^ 4^e^ c^^^ 



+1 "* o- , 1 i.<-^<^+i 



h' 2i + l r*^' 2i 

Integrating as before, we obtain for the potential of the 
vhole shell, 

(2j + 1)(2j + 3)/i* * 7'*^' 

Suppose now that we wish to obtain the potential of the 
whole shell at a point forming a part of its mass, distant r 
from the centre. We shall obtain this by considering sepa- 
rately the two shells into which it may be divided, the 
exteraal radius of the one, and the internal radius of the 
other, being each r. Writing r for a, in the first of the fore- 
going results, we obtain 

2i + l h' ^"^ ''^'^' 
And writing r for a in the other result, we obtain 

47r pP^T^^^-a!''^' 

(2i + l)(2i + 3) K 7*^^ • 

Adding these, we get for the potential of the whole 
sphere 

2i-hi A' I 2 ^■^(2i+3)r*^Y 
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It is hardly necessary to observe that the corresponding 
results for a solid sphere may be obtained from the foregoing, 
by putting a' = 0. 

If the density, instead of being ^^<c', be ^P^dl^y similar 

reasoning will give us, for the potential of the thin shell of 
radius c and thickness dc at an internal and external point 
respectively, 



And, integrating as before, we obtain for the potential of 
the whole shell, 

47r p a"*"^*"^' — a'*"""*' 
■(2i + l)(m + i + 3)r^* ?** ** ^ ^^*''°*^ F°^*- 



And, at a point forming a part of the mass, 
Si + l'FA m-i + 2 ^' "^ m-fi + 3 r 



11. Suppose, for example, that we wish to determine, in 
each of the three cases, the potential of a spherical shell 
whose external and internal radii are a, a', respectively, and 
whose density varies as the square of the distance from a 
diametral plane. 

Taking this plane as that of xy, the density may be ex- 
pressed by Ti^^, or pc^. Now jj? = — ^ — , Hence the 
density of this sphere may be expressed as 

The several potentials due to the former term will be, 

2 

writing 2 for i and multiplying by -j: , 
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1.) h* ^* ""^ ''^' 105 A«^« 7^ ' 15 h' \ 2 ^ 7r' J' 

And for the latter term, writing for i, and 2 for m, and 

inultiplying by ^ , 

t> 

*^^r.* n^^ ^iL?lr«'" 47r p / g*-/ r^-a^ 
j^-^^(,a-a;, ^^ ^, r"'^ 3 A«V 4 "*" 5r"r 

And, since P^r^ = — ^ — > ^® S®^ ^'^^ *^^ potential at an 
ntemal point 

it an external point 

at a point forming a part of the mass 

12. We may now prove that by means of an infinite series 
of zonal harmonics we may express any function of fi what- 
ever, even a discontiDUous function. Suppose, for instance, 
that we wish to express a function which shall be equal to 
A from /A = l to /a = X, and to B from /i=\ to ft= — 1. 
Consider what will be the potential of a spherical shell, 
radius c, of uniform thickness, whose density is equal to A 
for the part corresponding to values of fi between 1 and X, 
and to B for the part corresponding to values of [i between \ 
and « 1. 

Divide the shell, as before, into indefinitely narrow strips 
l>y parallel planes, the distance between any two successive 
planes 'being cdfi. 
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We have then, for the potential of such a sphere at 
point of the axis, distant z from the centre, 

for the first part of the sphere 






X y 



and for the latter part 

fzc r —4ft 



27r5c'^- ' '^^ 



These are respectively equal to 
27r^c''Sc 



:£ ^P. +P, ^ + P, J + ... + P,^l+ ....)c?/t, 



27r5c*Sc f* 



j\{Po + P^l + P>^+^" + P4+ -)<^l^. 



at an internal point ; and to 
27rAc'Sc 



2TrB(?hc t^ 



£(P, + P,^+...+P,^,+ ....)d/., 



at an external point. 

Now it follows from Chap. ii. (Art. 23) that if i be a 
positive integer, 

£p.<?/* = -^j{P^,(X)-P,.(X)l. 

whence, since I P^dll = 0, it follows that 
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^lao j P^dfl = l-\ j F,d/JL=l+\. 

i^nce the above expressions severally become : 
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the potential at an internal point on the axis 

c 

.2 



-^--^{P,(X)-P.(X)}^-.. 
for the potential at an external point on the axis 



z 
A-B 



{P.(X)-P.(X)}^-... 



Hence the potentials at a point situated anywhere are 
•ectively 

-^-f^{i'.W-P.(x)}^^ 

-^{P3(X)-P,(X)}^^-.... 

ia internal point; 
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aud 



Sttc'Sc [{A (1 - X) + B(l + X)} 



P.(/*) 



A-B 



3 ?{P,(X)-P.(X)}^^ 



at an external point. 

Now, if we inquire what will be the potential for the 
following distribution of density, 

^[4(1 -X) +B(1 + X)- (^-5){P,(X) -P„(X)}P» 

-(4-5){P,(X)-P.(X)}P.(/t)-... 
- {A-B){P^,i\) -F,.,{X)}P,Oi) - ...], 

we see by Art. 6 that it will be exactly the same, both at 
an internal and for an external point, as that above in- 
vestigated for the shell made up of two parts, whose densities 
are A and B respectively. 

But it is known that there is one, and only one, dis- 
tribution of attracting matter over a given surface, which 
will produce a specified potential at every point, both ex- 
ternal and internal. Hence the above expression must 
represent exactly the same distribution of density. That is, 
writing the above series in a slightly ditFerent form, the 
expression 

^ - ^ [^ + {P.(S) - PMP.M 

+ {P.(X)-P,(X)}P» + ... 
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• 

w equal to -4, for all values of /a from 1 to \ and to B for all 
values of /tt from X to - 1. 

13. By a similar process, any other discontinuous function, 

J^hose values aire given for all values of fi from 1 to — 1, may 

oe e^ressed. Suppose, for instance, we wish to express a 

Auction which is equal to A from fi = l to fi=^\, to n from 

A**Xj to /i=sX^ and to Ofroixi /a =X^ to /a=^ — 1. This will 

be obtained by adding the two series 

For the former is equal to -4 — J5 from /i = l to /A = Xp 
«id to from /i = X^ to /^ = — 1 ; and the latter is equal to 
B.from /i6 = 1 to /i6 = Xj, and to C from /*.= X, .to jm =— 1. 

By supposing A and (7 each = 0, and B = 1, we deduce a 
Beries which is equal to 1 for all values of ft from iA = \to 
M=Xy and zero for all other values. This will be 



This may be verified by direct investigation of the 
potential of the portion of a homogeneous spherical shell, 
of density unity, comprised between two parallel planes, 
listant respectively c\ and c\ from the centre of the 
'Pherical shell. 

14. In the case in which X^ and X^ are indefinitely nearly 
qual to each other, let X, = X, and Xj = X + dX, We then 
^ve, ultimately, 



i^A)-^AJ=^^ 
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Hence P^,(\) - P^,(X^) -P^»(\) - P^.(\) 

= (2i+l)P,(X,)(BL 
Hence the series 

^ {1 + 3P;(X)P» + 5P,(X)P» + ... 

+ (2i + 1) P,(X)P,(a*) + 

is equal to 1 when /Lt = X (or, more strictly j when fi has 
value from X to X + dX) and is equal to for all other vali 

V\ 5^5 s c. - We hence infer that 

'^ S i^ 1+ 3P,(X)P» + ... + (2i + l)P,(X)P,Oi) + ... 

^ H ( is infinite when /a = X, and zero for all other values of /Lt. 

15. Bepresenting the series 

. \[1 + 3P,(X)P,(/i) + ... + (2i+ 1)P,(X)P,0*) + ...} 

by ^(\) for the moment, we see that p^(X)c2X is equal to ^ 
when fi^\ and to zero for all other values. Hence tkv 
expression 

{Pi<^(\)+P«<^(\) + -}^ 

is equal to p^ when fi = X^, to p, when /it = >c,... Supposmg 
now that X^, X^... are a series of values varying continuous^ 
from 1 to — 1, we see that this expression becomes 

/ p<^(X)dX, 

p being any function of X, continuous or discontinuon& 
Hence, writing <^(X) at length, we see that 

||£ peZX + 3P»|' pP,(X)(ZX+... 

+ (2i + l)P,(/*) J* pP,(X)(iX + ...| 

is equal, for all values of fi from — 1 to + 1, to the Bame 
function of fi that p is of X. 



TO THE THEOBY OF ATTRACTION. 67 

16. The same conclusion may be arrived at as follows : 

The potential of a spherical shell, whose density is p, 
kiid volume U^ at any point on the a:ds of z, is 

pdK 






which is equal to g-j I pd\+-j pPi(\)cZ\+ ... 

for an internal point, 

and to ^i^f pdK'{-^J%P^{\)d\+... 

^^j j)P[{\)d\ + .X 

for an external point. 

It hence follows that the potential, at a point situated 
anywhere, is 

^[jd\ + ?^^jjP^iX)dK + ... 

for an internal point, 

tad to ji^f pd\ + ^^f pP.OOdk + . . . 

or an external point. 

And these expressions are respectively equal to those 
or the potentials, at an internal and external point re- 
pectively) for matter distributed according to the following 
awof den£iit/: 

lb— ^ 
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+ {2i+l)P,Qi)f pP,(K)d\-\- J. 

It will be observed, in applying this formula, that if pbe 
a discontinuous function of \, each of the expressions of the 

form I pP^{\)dK will be the sum of the results of a seriies (tf 

integrations, each integration being taken through a series of 
values of X, for which p varies continuously* 



CHAPTER IV. 

^HERICAL HAKMONICS IN GENEEAL. TESSERAL AND SEC- 
TOBIAL HARMONICS. ZONAL HARMONICS WITH THEIR 
AXIS IN ANY POSITION. POTENTIAL OF A SOLID NEARLY 
SPHERICAL IN FORM. 

1. We have hitherto discussed those solutions of the 
quation V'F= which are symmetrical about the axis of 2?, 
r in other words, those solutions of the equivalent equation in 
olax co-ordinates which are independent of ip. We propose, 
X the present Chapter, to consider the forms of spherical 
annonics in general, understanding by a Solid Spherical 
[annonic of the i^ degree a rational integral homogeneous 
motion of X, y, z, of the i^ degree which satisfies the equa- 
oa V*F= 0, and by a Surface Spherical Harmonic of the 
* degree the quotient obtained by dividing a Solid Sphe- 

^ Harmonic by (55* +y* + «*)*• Sudb an expression, as we 
e by writing «=rsin^cos^, y = rsin^sin^, z^rcosd^ 
H be of the i*^ degree in sin cos ^, sin sin <f>, cos ; and 
11 satisfy the diflferential equation in Y^ 

Id/. .dY,\^ 1 cPF, .,.^-.^ ^ 
► writing fM for cos 0^ 



U'-'^ 



f}+r^-^+'('+"J'.-«- 



It will be convenient, before proceeding to investigate the 
gebraical forms of these expressions, to discuss some of 
leir simpler physical properties. 

2. We will then proceed to shew how spherical har- 
onics znajr be employed to detenmxie \5dl^ ^X^oSy^^ %sA 
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consequently the attraction^ of a spherical shell of indefinitelj 
small thickness. 

We will first establish an important theorem, connecting 
the potential of such a shell on an external point with that 
on a corresponding internal point. The theorem is as follows; 

I f O he the centre of such a shelly c its radius^ P any «n- 
temal point, F an external point, so situated that Y lies on 
OP ^^oduced, and that OP . OF = c', and if OP = r, OF =1^, 
th/en the potential of the shell at P is to its potential at F 
asctoT^or {which is the same thing) asi' to c 

For, let A be the point where OP* meets the surfeoe of 
the sphere, Q any other point of its surface. Then, by a 
known geometrical theorem, 

QP: QP v.APiAF v.C'-ri /-a 

. , c — r cr — ♦^ cr-'f^ r o 

And -7 — = — ? = ij = -«=:-,. 

r — rr ^ cr o —cr o r 

Again, considering the element of the shell in the im- 
mediate neighbourhood of Q, its potential at P is to its 
potential at P as QP is to QP, that is, as c to r, or (which 
IS the same thing) as r' to c, which ratio, being independent 
of the position of Q, must be true for every element of Ihe 
spherical shell, and therefore for the whole shelL Hence 
the proposition is proved. 

3. Now, suppose the law of density of the shell to he 
such that its potential at any internal point is FQi, if) -n* 

r* 

Then F{/Jb, ^) -^ must be a solid harmonic of the d^:eei 

c 

Hence FQi, ifi) must be a surface harmonic of the d^ree u 

Let us represent it by Y^. 

By the proposition just proved, the potential at any 
external point, distant r from the centre, must be 

Y ~— 
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Hence, the component of the attraction of the sphere on 
3 internal point measured in the direction from the point 
WQxdB, ie. towards the centre of the sphere, is 

And the component in the same direction of the attraction 
I the external point, measured inwards, is 

Now suppose the two points to lie on the same line 
issing through the centre of the sphere, and to be both 
definitely close to the surface of the sphere, so that r and r 
6 each indefinitely nearly equal to c 

And the attraction on the external point exceeds the 
•traction on the internal point by 

(2t + l)-5. 

Now, supposing the shell to be divided into two parts, 
J a plaiie passing through the internal point perpendicular 
) the line joining it with the centre, we see that the at- 
•action of the larger part of the shell on the two points will 
B ultimately the same, while the component attractions of 
le smaller portions, in the direction above considered, will 
3 equal in magnitude and opposite in direction. Hence the 

fference between these components, viz. (2i+l) — , will be 

c 

[Ual to twice- the component attraction of the smaller 

•rtion in the direction of the line joining the two points. 

It if P( be the density of the shell. Be its thickness, this 

inponent attraction is 2irpiSc. 

Y 
Hence (2* + 1) — * = 4t7^p^ 8c, 

c 

_.2t + l ^ 
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And, if a^ be the corresponding surface density, 

2i + l 



€r,=: 



r. 



It hence follows that if the potential of a sphericalsheB, 
of indefinitely small thicmees, be a surface harmonic, its 
potential at any internal point will he proportional to ih$ 
corresponding solid hamionic of positive degree, and its po- 
tential at am,y eoctemal point will he proportional to the 
corresponding solid harmonic of negative degree. 

That is, the proposition proved fot zonal harmonics in 
Chap. m. Art. 6, is now extended to spherical harmonics in . 
general. 

I 

4. The spherical harmonic of the degree i wiU irwclve 
2i + 1 arbitrary constants. 

For the solid spherical harmonic. r'Y], being- a rational 
integral function of x, y, z of the ir degree, will consist of 

(i + l)(i + 2) ^ TVT XI. ' • T^lTr V • 

Y ' terms. Now the expression V F, bemg a 

rational integral function of x, y, z of the degree t — 2, will 

consist of "^ — 5-^ terms ; and the condition that it must be 

(i—VSi 
= for all values of x, y, z, will give rise to ^ — ^ — relations 

among the ^ ^-~ coefficients of these terms, leaving 

^^ -^ ' — - — ^^ , or 2i + 1, independent coefficients. 

6. We proCfeed to shew how the spherical harmonic of the 

degree i may be arranged in a series of terms, each of which 

may be deduced by diflferentiation from the Zonal Harmonic 

sjnnmetrical about the axis of z. The solid zonal harmonic, 

which, in accordance with the notation already employed, is 

represented by r^P^ (jjl), is a function of z and r of the degree », 

d^V d^V d^V 
satisfying the equation V* F= 0, or -7-^ + -j-j + -tt = 0- 

Now, i{ we denote ttiis ^•sspces^oTi \j^ P^V^^^^^^ftathat 
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ice it is a function of z and r, it is a function of the dis- 

uice («) from a certain plane passing through the origin, and 

E the distance (r) from the origin. Further, if we write for z 

be distance from any other plane passing through the origin, 

, , . d^'V cPV dPV ^ .„ 

eaving r imaltered, the Equation -r-j +;j^ + '^ =^ '^"l 

30iitmue to be satisfied. 

Now z-Vaix+J-^ly^ a being any quantity whatever, 
represents t^e distance f^om a certain plane passing through 
the origin, since in this tepressioA, the sum of the squares 
)f the coefficients of z, x, y is equal to imity. Hence 

Pi(« + a(a? + >/— ly)} is a solid zonal harmonic of the 

on 'U 

legree t, its axis being the imaginary line - = /= = ^• 
lierefore the equation 



satisfied by F= P, {« 4- a (a? + V— 1 y)}, that is, expanding 
f Taylor's Theorem, it is satisfied by 

/N. / . /— T sdPA^y a' ^ i—= .^d^P,{z) ^ 

of(x^^r:iyYdT,{z) 

"^ 1.2... i dz' ' 

^ aU values of ou 

Hence, since the equation in V is linear, it follows that 
is satisfied by each term separately, or that, besides P^ {z) 
self, each of the i expressions, 

.wrx»'?|!),(«+,/ziy).^....(«+V^,).*p, 

^tigfies the equation F=0. 

By similar reasoning we may shew that each of the i ex- 
^688ions, 

tiafies iihe same equation. 
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Now each of the 2i solutions^ thus obtained, is imaginaiy. 
But the sum of any tw6 or more of them, or the result 
obtained by multiplying any two or more by any arbitraiy 
quantities, and adding the results together, will also be a 
solution of the equation. Hence, adcung each term of the 
j&rst series to the corresponding term of the second, we ob- 
tain a series of i real solutions of the equation. Another 
such series may be obtained by subtracting each term of the 
second series £rom the corresponding term of the first, and 

dividing by V— 1. We have thus obtained (including the 
original term Pi{z)) a series of 2i*+l independent solutions 
of the given equation, which will be the 2* + 1 independent 
solid harmonics of the degree t. 

6. We may deduce the surface harmonics from these by 
writing r sin 5 cos ^ for a?, r sin sin ^ for y, r cos for «, 
and dividing by r*. Then, putting cos = fjb, and observing 

that P,(«) = r*P,0*), ^^)-:=r*^i^... we obtain the fol- 

lowing series of 2i+ 1 solutions : 

cos^sind^^, cos2^sin»^'^^^, ...coai<l>Bm*0^^, 

Expressions of the form 

Ocosff^sin'^^^l^, 

or asm. 0-9 sin'0 ,^ , 

or their equivalents, 



(7C08^^(1-Mf.^^\ 
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ind 8 denoting any quantities independent of 6 and <}>) 
called Tesseral Surface Harmonics of the degree i and 
r a. The particular forms assumed by them when 
i are called sectorial Surface Harmonics of the degree i. 

ill be observed that, since J \* is a numerical constant, 

orial Harmonics only involve 6 in the form 

sin*^,or(l-/x,«)»\ 

Phe product obtained by multiplying a Tesseral or 
orial Surface Harmonic of the degree % by r* (that is, 
expression directly obtained in Art. 5) is called a Tesseral 
3ctorial Solid Harmonic of the degree i. 

\ We shall denote the factor of a Tesseral or Sectorial 

d'P (a) 
nonic which is a function of 0. that is sin*^^ — j^^^ , or 

- d'P (a) 
/Lt*)* js 9 * ^y ^® symbol 7J(*^>, or, when it is necessary 

articularize the quantity of which it is a function, by 
(/i) or T/') (cos 6). 

t will be convenient, for the purpose of comparison with 
forms of Tesseral Harmonics given in the Micanique 
tie, and elsewhere, to obtain Tf"^^ in a completely de- 
)ed form. 

low, since P^Qi) - q* i 2 3 — * — d *" ' ^® ®^ *^* 

[fjL)_ 1 d*^^(fi'--lY 

' "■2M.2.3...* rf/A»+«' 
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- J (2i - 2) (2i - 3)...(i - <r - 1) /t<— « 

+ ^-^Y~^ (^*~'*^ (2»-5)...(i-<r-3)At*-'-* 



(^-a^)(^-«^-I)(^-«^-2)(t-;^-g) ^,^ | 
■^ 2.4.(2i-l)(2t-3) '^ "'T 



And therefore 



2t(2t-l)...(w+l) I f ._. (wXi-^1) ^:... 

^' - 2M.2.3...* ^^ '*^ Y a(2»-l) '* 



(i_o-)(t-«r-l)(*-o--2)(«-o-3) 
+ 



2.4(2i-l)(2t-3) 



\ u.*-'-*— I 



The form given by Laplace for a Tesseral Surfeu^ Hitf>- 
monic of the degree i and order cr is (see MSoanique CHede, 
Liv. 3, Chap. 2, pp. 40—47) 

A being a quantity independent of 6 and ^ The factor of 
this, involving fi, is denoted by Thomson and Tait (Ifatuni 
Philosophy, Vol. 1, p. 149) by the symbol ©/'>. Thomaon 
and Tait also employ a symbol ^/^\ adopted by Maxwell in 
his Treatise on Electricity and Magnetism, Vol. 1, p. 164 
which is equal to 
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Seine represents the expression 

(t-- <r) (t--«r - 1) (» - «r - 2) (t - «r - 3) ) 

"^ 2.4.(2t-l){2»-3) '^ •••)•■ 

[— 1)* ©/*^>, by the symbol PJ(jj), and calls these expres- 
18 by the name Zv^eordnete Functionen Master Art {Hand- 
h der Kugelfumtionen, pp. 117, 118) which Todhunter 
cislates by the term ''Associated Functions of the First 
id,'* which we shall adopt. 

Heine also represents the series 

^ (t-(r)(i-<r-l) ,,,,, 
2(2<-l) ^ 

■*" 2-4(2i-l)(2i-3) '^ 

the symbol ^^{/i), (p. 117). 

The several expressions, 1Y\ ©1'), ^'>, PI, ^U are con- 
:;ted together as follows : 



(2i -!)...(» -«r+l) 



yw = (») 



_ 2*"**(»-l)...(g- + l) cv/rt _ /_ i\i pf-d- uV 39' 
''(t + <r + l)(i+o- + 2)...2t^ ^ ^^ ^' ^^ '^^ *" 

8. It has been already remarked that the roots of the 
raation P. = are aJl real. It follows also that those of the 

dP d^P 

luations ^ = 0, -r-j ^ 0... Are real also. Hence we may 

"rive at the following conclusions, concerning the curves, 
aced on a sphere, which result from our putting any one 
these series of spherical harmonics = 0. 

By putting a zonal harmonic =0, we obtain t small circles, 
lose jAanea are parallel to one anotlieT, -jerjoxi^vsciNax \ft 
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the axis of the zonal harmonic^ and symmetrically situated 
-with respect to the diametral plane, perpendicular to this 
axis. If i be an odd number this diametr^ plane itself 
becomes one of the series. 

By putting the tesseral harmonic of the order o-^O, we 
obtain i — a- small circles, situated as before, and a great 
circles, determined by the equation coscr^ = 0, or siii<r^ = 0, 
as the case may be, their plaiies all intersecting in the axis 
of the system of harmonics, the angle between the planes of 

any two consecutive great circles being - . 



By putting the sectorial harmonic =0, we obtain i 
great circles, whose planes all intersect in the axis of the 
system, the angle between any two consecutive planes being 



TT 



9. The tesseral harmonic may be regarded from anotha 
point of view. Suppose it is required to determine a soUd 
hannonic of the degree t, and of the form F^r*, such that F, 
shall be the product of a function of fi, and of a function of ^ 
which functions we wUl denote by the symbols Jlf^, <E>f, respec- 
tively. The differential equation, .to which this will lead, is 

Now this will be satisfied, if we make M^ and <E>| satisfy 
the following two equations : 

The latter equation gives 

O, = C7cos <r(f> + C sin <r^. 

And, faking cr as an mle^et, ^«v!w«^ ot ue^tire, the 
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former is satisfied hj M,^ T^(^\ le. {l^p^ifj-j (l-^fij, 
as we proceed to prove. 

We know that 

Differentiate o- times, and we get 
whence, by Leibnitz's Theorem, 

t- or 

andj multiplying by (1 —/a')*, 

, +(.--<r)(t+,r+l)(l-/*V0'=O...(l). 

Now, putting av/^'^y/". 
we get 

£p^» (1 - ^^f J!;!? _ .^ (1 - ^«)-'-' ^. . 



-f( 
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And «(» + !) r/-) = i(,-+X)(l-;r^^|5. 

I 

+ (t-<r) (t + <r + l) (l-/*')^P'+a'(l -/*«)^-*^ 

=«r«(i-/*T^*by(i) 

1 — /I. 

Hence the equation above given for 3f^ is satisf 
if, = jT/'), and the equation in F, is satisfied by 

y; = Cr/') cos (T^ + G' T/') sin cr^ 

10. In Chap. Ii. Art.lO we have established the funda 
property of Zonal Harmonics, that if i and m be two u 

positive integers, I P^P^dfi = 0. This is a particul 

of the general theorem that if TJ, 5^ be two sar£eM 
monies of the degrees i and m respectively, 



|'J^*'r.r„(£,«z<^=o. 
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^% let T^ F^ be the corresponding solid harmonics, so 
wiat V^=.r*Y^, 7^ = 7^ Y^. Then, by the fundamental pro- 
i^rtjr of potential functions, we have at every point at which 
^ attracting matter is situated, 

da?'^ djt'^ M ' da?~^ df "*" da* ""' 
>nd therefore 

r, in accordance with our notation, l^v'^"" 1^V*^ = 0. 

Now, integrate this expression throughout the whole 
pace comprised within a sphere whose centre is the origin 
nd radius a, a being so chosen that this sphere contains no 
ttracting matter. We then have 

l\\{y,v'V^-V^^'V:idxdydz = 0. 

But also, when the integration extends over all space 
Jomprised within any closed surface, we have 

\\\{V;^V^- F^v'T^ dxdydz=jj(r, 5 - F„^') d8=0, 

is denoting an element of the bounding surface, and -r- 
lifferentiation in the direction of the normal at any point. 

Now, in the present case, the bounding surface being a 
tphere of radius a, and FJ, V^ homogeneous functions of the 
legrees i, m, respectively, 

nd, the integration being extended all over the surface of 
be sphere, the limits of /4 are — 1 and 1, those of <^, and 27r. 
lence 
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whence, if m--i be not ^0, 



u. 







/I r2ir 
I Y^dfid(l> will be investigated here- 
after. 

11. We may hence prove that if a function of ii wnA ^ 
can he developed in a series of surfdce harmonics, ^uch (fc- 
veloprnent is possible in only one way. 

For suppose, if possible, that there are two such develop- 
ments, so that 

FQi, ^) = F,+ r; + ... -f F, -f ... 

and also 

Then subtracting, we have 
0= r^- F;+ F,- F/+ ... + F,- T; + ... identically. 

Now, each of the expressions F^ — F/, i^ — F/... F,- 7/ 
being the diflference of two surface harmonics of the degree 
0, 1, ...i... is itself B, surface harmonic of the degree 
0, 1, ...*.... Denote these expressions for shortness by 
Z^, Z^.^.Z^ ..• so that 

0=-^ + -^ + .,. +-2] + ... identically. 

Then, multiplying by Z^ and integrating all over the 
surface of the sphere, we have 



0=1 f' Z^^dfidij). 



That is, the sum of an infinite number of essentially 
positive quantities is = 0. This can only take place when 
each of the quantities is separately = 0. Hence Z^ is identi- 
cally = 0, or F/ = Y^, and therefore the two developmentB 
are identical. 

We have not assumed here that such a development is 
always possible. That it ia so, Nr!S\.\i^ ^^'wa.V^\«o&ftx^ 
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12. By referring to the expression for a surface har- 
monic given in Art. 4,'i7e see that each of the Tesseral and 
Sectorial Harmonics involves (1 — /**)*> or some power of 

(1 - ij!)\ as a factor, and therefore is equal to when /* = + !. 
From this it follows that when /* = ± 1, the value of the 
Surface Harmonic is independent of </>, or that if T (ji, <l>) repre- 
sent a general surface harmonic, Y(±l, if>) is independent of 
^ and may therefore be written as F (+ 1). Or Y (1) is the 
Value of YQi, ^) at the pole of the zonal harmonic P, (/i), 
F(— 1) at the other extremity of the axis of .P, (ji)^ 

We may now prove that 



L 



r^-2,rr.(i)i>,M. 



For, recurring to the fundamental equation, 

Now, if we integrate this equation with respect to <^, 
etween the limits and 27r, we see that, since 



/ 



d*Y dY, 



nd the value of Y^ only involves if> under the form of cosines 
r sines of d> and its multiples, and therefore the values of 

TT are the same at both limits, it follows that 



Jo W^"^"^^' 



Hence 



: |i{a-/*')(/rw)}+»-(i+i)(f>.#)=o. 

Hence I rj^<^ is a function of fi which satisfies the 
Jo 
iindamental equation for a zonal harmonic, and we therefore 

ave 
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•a* 



r 

^0 





G being a constant, as yet unknown. 

To determine G, put /a=1, then by the remark just made, 
Y^ becomes Y^{l)y and is independent of if>. Hence, when 

r2ir 

/A = l, J r,d^=:27rr,(l). Also P,(/a) = 1. Wehaveihere- 
fore ^ 27rr,(l) = (7, 

/. rY,d^^2irY,{l)PM. 
Jo 

It follows from this that 

13. We may now enquire what will be the value of 

/I r2«- 
/ Y,Z,dfid<l>, 

F,, Z^ being two general surface harmonics of the degree i 
Suppose each to be arranged in a series consisting of the 
zonal harmonic P^ whose axis is the axis of z, and the system 
of tesseral and sectorial harmonics deduced from it. Let us 
represent them as follows : 

r,= AP, 

+ C;T/i> cos ^ + G^T^'^) cos 2^ + ... + C^T^^^^ cos aji + ... 

+ G.T^^emi^ 

+ 8JP sin + 8Ji^^ sin 20 + ... + S^TI^'^ sin o-<^ + ... 

+ >8;r,»sin«^; 
^,= aP, 

+ Cj T,(i> cos + Cj r/2) cos 20 + ... + c<,2;<^> cos ^0 + . . . 

Ca^i^*^ cos i^ 
+ 5,r,<i>sin0 + 5,y/2)sin20 + ...+5ay/*')sin^0 + ... 

+ «^r/*'>8ini;^. 

Hence the product Y^Z^ will consist of a series of terms, 

in which ff> will enter under the form cos 0-0 cos <r'0, or 

cos cr^ sin a'^. This expieasvoii ^\vfivi ydX*^^^^ laetween 
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the limits and 27r vanishes in all cases, except when 
<T' = <r and the expression consequently becomes equal to 
cosV<^, or sin\o-^. In these cases we know that, cr being any- 
positive integer, 

/ftr rair ^ 

cos' <7^ (2^=9 I SID? a(j>d^ = TT. 

Hence the question is reduced to the determination of the 
value of 



/: 



-1 



Now !r/') = a-/*? ^ 

But, by the theorem of Bodrigues, proved in Chap, II. 
&it. 8, we know that 

'^^ \iEz d/**-' • 

Hence 7/") may also be expressed under the form 

'^ ^ 2'.1.2.3.'..i[i-£'' '^'' d/*<-' ' 

rhence it follows that 



Now, putting (/t* — l)* = Jf for the moment, and inte- 
rating by parts, 
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The factor -r-; — vanishes at twth limits, hence 

by a repetition of the same process. 

And by repeating this process o- times, we see that 

].^-W' ■di^F^^'^^^-^^'j.^W) ^^ 

= (- 1)' (2M . 2 . 3 ... i7 P P/ d/* 

c 

= (-l)'(2M.2.3...t)*2^. 
and therefore 

It will be observed that this result does not hold whe 
o- =5 0, in' which case we have 

Hence j r'Y,Z,dftd<l> 



/t \i + <r 9 



47r J 



2i + l 



* In this case J ' toft* c4>d4>= r* ^x^ v^il^wix. 
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^ (li + l |«+2 



. 14, We have hitherto considered the Zonal Harmonic 
under its simplest form, that of a " Legendre's CoeflSicient " in 
^hich the axis of 0, ie. the line from which is measured, is 
the axis of the system. We shall now proceed to consider it 
luder the more general form of a "Laplace's Coefl5cient>' 
iu which the a^ of the system of zonal harmonics is in any 
position whatever, and shall shew how this general form may 
36 expressed in terms of P<(/a) and of the system of Tesseral 
ttid Sectorial Harmonics deduced from it. 

Suppose that ^, <^' are the angular co-ordinates of the 
ixis of the Zonal Harmonic, i.e. that the angle between this 
•xis and the axis of z is d^, and that the plane containing 
hese two axes is inclined to a fixed plane through the axis 
f z which we may consider as that of zx^at the angle <f)\ 
n accordance with the notation already employed, we shall 
epresent cos ^ by fi\ 

The rectangular equations of the axis of this system 
rill be 

X ^ y ^ ^ 

sin ff cos <^' sin ff sin <^' "" cos ^ ' 

Hence the Solid Zonal Harmonic of which this is the axis 
' deduced from the ordinary form of the solid zonal har- 
monic expressed as a function of z and r by writing, in place 
f «, a? sin ff cos ^' + y sin ff sin ^' + « cos ff. 

To deduce the Surface Zonal Harmonic, transform the solid 
>iial harmonic to polar co-ordinates, by writing r sin d cos ^ 
>^ a?, r sin 5 sin <^ for y, r cos for z, and divide by /. 

The transformation from the special to the general 
*lai of surface zonal harmonic may be at once eflfected, 
f substituting for /*, or cos 0y cos^cosa +sin^sin^cos(0— ^'). 

Now, in order to develope 

P^ [cos ^ cos ^ + sin ^ sia fl' cos (<^— <^"|\ 
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in the manner already pointed out, assume 

P^ {cos Ocosff + sin sin ^ cos {<f) — 0')} 

= ^P,0^) + (C(i)cos^ + iSfWsin<^) 2J1) 

+ (0(2)cos2<^+/S<^sin2^)ri2) + _^ 

+ ((7<'>coscr^+/8^<'')sin<7^)r/<^)+... 
+ (C(»^cosi0 + >8fWsini»27*^ 

the letters u4, ... C^'\ 8^^^... denoting functions of /fraud 
if), to be determined. 

To determine C^*^\ multiply both sides of this equation 
by cos a<f>T^^''^ and integrate all over the surface of the sphere, 
i.e. between the limits — 1 and 1 of ft, and and 27r of f 
We then get 

r [ 'P, {cos cos ff + sin sin ff cos {<f> - <^')} <5os <r^ 2J''> cW 

= C W r r (cos ^^j;<<'))* d^^i^ 



[i_+€P 27r 



|^-o• 2^4•l 



C(-). 



It remains to find the value of the left-hand member of 
this equation. 

Now cos o-^jy*^) is a surface harmonic of the degree i, and 
therefore a function of the kind denoted by T] in Art. 12. 

And we have shewn, in that Article, that 

that is, that if any surface harmonic of the degree i be muUi- 
plied by the zonal harmonic of the same degree, and the product 
integrated all over the surface of the sphere, the integral is 

eqtuil to ^. — Y into the value which the surface harmonic 

assumes at th^pole of the zonoX Karmomo. 



. 
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Mii therefore 

j J Pi {cos ^ cos ^ + sin ^ sin ff cos (<f> - <^')1 <^s <r4>Ti'^dfid<l> 

Hence 



or C<') = 2 14=^ cos ad> T^^^ (ji). 
\% + a' 

Similarly iSfW = 2 ^t^ sin erf 2;(<^) (/li')- 



i + <r 



And to determine A, we have 
^ {^P^ {cos^cos^ + sin^sin ^cos (^- 0')1 Pi W ^/^# 

or A = P,{fi')' 
Hence, P, {cos (? cos ^ + sin ^ sin ^ cos (^ — <j>)} 

-na*')P,0t)+2J=co8(^-f)7;a)0*')2'.«0*) 



+ 2 



Jyl COS 2 (0 - ^') r.® 0*') 2;® 0*) + 
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r * u \ 

+ 2L^cos<r(^-<^')2',W(/t')2'/'>0*) + ..- f 

1* + *^ ll 

+ 2 j^ cos i (^ - f ) r.® (^') r.w H 

15. We have already seen (Chap. n. Art. 20) how any 
rational integral function of fi can be expressed by a finite 
series of zonal harmonics. We shall now shew how any 
rational integral function of cos^, sin cos ^, sin ^ sin ^ 
can be expressed by a finite series of zonal, tesseral, and 
sectorial harmonics. 

For any power of cos ^ or sin ^, or any product of sucli 
powers^ may be expressed as the sum of a series of terms of 
the form cos o-^, or sin a^, the greatest value of c being the 
sum of the indices of cos^ and sin ^, and the other values 
diminishing by 2 in each successive term. Hence any 
rational integral function of cos 0, sin cos ^, sin sin <f>, will 
consist of a series of terms of the form 

cos*^ sin" cos a(f> or cos*" sin" sin a<j>, 

where n is not less than <r. 

If n be greater than cr, w — a- must be an even integer. Let 
71 — <r = 25, then writing sin"^ under the form (1 — cos* 5)' sin'tf, 
we reduce cos*^ sin" cos <7^ to the sum of a series of terms 
of the form cos^ &m^ cos a^, or, writing cos 0^ fi,oi the 

form fi^ (1 — /**) * cos aij>. 

Similarly cos*" ^ sin" ^ sin <70 is reduced to a series of 

terms of the form /jF (1 — /x") * sin o-^. 

and /M^"*"' can be developed in a series of terms of the form 
of multiples of Pp+^, JPp+^-2 •.. . (Chap. n. Art. 17.) 

Hence fiP can be expressed in a series of the form 
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A^, A^ representing known numerical constants^ and therefore 
M*^ (1 — fi^y assumes the form 

consequently multiplying these series by cos o-^ or sin o-<^, we 
obtain the developments of 

/aP(1 — /x")* coso-^ a-nd /aP(1 — /x')*sino-0 
in series of tesseral harmonics. 

16. We will give two illustrations of this transformation. 

First, suppose it is required to express cos* 6 sin' d sin <^ cos ^ 
in a series of Spherical Harmonics. 

Here we have sin ^'cos = ^ sin 20. 

Hence cos* sin* sin cos ^ = ^ cos* sin* sin 2^. 

Comparing this with cos" siu" sin <r<f>, we see that n is 
tot greater than a. 

Hence cos* d sin* sin <]>eoa<f> = ^(j?{l — /**) sin 2<t>. 

•* '^ ~12V35 dn' "^ 7 d/i' J 

105 d/ "*" 21 d/t"' 
cos* ^ sin* d sin ^ cos tf> 

= lJi_€§(l-^») sin 2^ + 10(1-/**) sin 2^} 



2 1105 dfi" 
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Next, let it be required to transfonn cos*^ sin* sin ^ cos*^ 
into a series of Spheric5al Harmonics, 

Here sin <j> cos^ ^ = ^ sin 20 cos = t (sin 30 + sin <j>). 
Now cos' sin» ^ sin 30 = /*' (1 - /a*)^ sin 30 

Also cos' sin' ^ sin ^ = /it* (1 — /*') (1 — /t*)i sin ^ 

= (At'-/it»)(l-/it')isin^ 



Also (Chap. n. Art. 17) 

^•-16 p . 24p . lOp .Ip 
'^ "231 •"*'77 * 21 » ? •* 

Hence cos' sin* ^ sin 3^ 

= 120(23-1^ +77^]^^-''^'''^^* 

= {3^5^.''^ + 34^^-3^- 
Andcos'^sxn'^sm^ = -(gg3^' + 7^^* + g3-^ 

/ 8 rfP, 2 dP^ 4 dP,\ ,_ ^1 . . 
= -(693-^-385 ■^-B3^j(^-'*^*^* 

= - T-— rW- — TO)- — T'^A sin 6 : 
V693 * 385 * 63 * )^^9f 

.-. cos'^sin*^ sin^ cos*0 = |g^ r.W + j^ rwJsinS^ 

_ I-?_ TW _ JL r a) — L TO)! sm4. 
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!?• The process above investigated is probably the most 
mvenient one when the object is to transform any finite 
igebraical function of cos 6, sin cos <^, and sin 6 sin ^, into 
series of spherical harmonics. For general forms t)f a 
motion of ^ and ^, however, this method is inapplicable, 
ad we proceed to investigate a process which will apply 
niversally, even if the function to be transformed be discon- 
nuous. 

We must first discuss the following problem. 

To determine the potential of a spherical shell whose 
irface density is F(jx, ^), JP' denoting any function whatever 
f finite magnitude, at an external or internal point. 

Let c be the radius of the sphere, / the distance of the 
oint from its centre, ff^ (f>' its angular co-ordinates, V the 
otential. Then fi being equal to cos 

J-iJo [r'^-2cr'{cos^ cose'+sin^sine'cos(^-0')} + c']* ' 

The denominator, when expanded in a series of general 
onal harmonics, or Laplace's coefficients, becomes 

>T an internal and an external point respectively, P, (/*, ^) 
eing written for 

P, {cos ^ cos ^ + sin sin 0' cos {(f) — 0')}. 

Hence, F^ denoting the potential at an internal, F, at an 
eternal, point| 

; ^oifJ'^FOi, <l>)d,id<l> + ^j'^J'^P,{fL, 0) F(jM, 4>) dfid<f> 
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It will be observed that the expression Pi{fi, 4>) ^^^^ 
fA and fi symmetrically, and also <^ and ^\ Hence it sati 
the equation 

And, since fi and 6 are independent of fi and j!, 
differential equation will continue to be satisfied after jP|! 
been multiplied by any function of fA and <f>, and integn 
with respect to fi and ^. That is^ every expression of i 
form 

is a Spherical Surface Harmonic, or "Laplace's Functiott 
with respect to fi and ^' of the degree i. And the sevi 
terms of the developments of V^ are solid harmonics of 
degree 0, 1, 2...i... while those of F, are the correspo: 
functions of the degrees —1, —2, — 3.., — (i+ 1), ... — . 
these are the expressions for the potential at a point (r\ /4,flj 
of the distribution of density F(jm\ ^') at a point (c, /*', ^J. 

Now, the expressions for the potentials, both extendi 
and internal, given in the last Article, are precisely the saitt' 
as those for the distribution of matter whose surface densityii 

+ (2i + 1) ff'PiO^> 0) F(fi, <l>) dfid<f> + .. j, 
or, as it may now be better expressed, 



ilfS^^-*'^^'^ 
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3 J j P^ {cos^ cos^'+ sin^ sin^ cos (^-^') F(jiy if>) dfjud^ 






h 



And, since there is only one distribution of density which 
nil produce a given potential at every point both external 
knd internal, it follows that this series must be identical 
nth Fiji', 40' We have thus, therefore, investigated the 
tevelopment of F(ji', <l>) in a series of spherical surface 
utnnonic9*. 

The only limitation on the generality of the function 
^(m > <^') is that it should not become infinite for any pair of 
ulue$ comprised between the, limits —1 and 1 of fi, and 
knd 27r of ^. 

18. Ex. To express cos 2^' in a series of spherical har- 
nonics. 

For this purpose, it is necessary to determine the value of 
2i + 1) J f V, {cos^ cos^+sin5 sin^ cos (0-^')} c^s 2^dfid<f). 

Now Pt {cos ^ cos ^ + sin sin ff cos (<j> — 0')} 

= P,(cos5)P,(cos5') 

^ 2 . ^dP,(cos0) . ^dP, (cos ff) .. .,, 

€?P,(C0S^) . a^(?P,(C0S^) o/. .n , 

' df.^ ^ "^^ d/.-^ ^ C0S2(<^^^) + ... 

Now / COS €r (^ — <f>) cos 2^ d!^ = 0, 

or all values of a except 2. 

* In 'eonnection with the sabjeot of this Article, see a papor by Mr G. H. 
Hurwin in the Messenger of MathenuLtia for Maxoih) 1^11 • 
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And I cos2(^ — ^'}cos2^(2^ssTrcoB2^. 
Jo 

Also 
And 

ss2M • 2 • S.-.z/x^Pi g-^=i • 

Now when /x = 1, 

And when /it=» — 1, 

Hence 

£sin'^g^^ = ^,^j3 . 2«.1.2.3...qi-(-in 

= 4 or 0, as i is even or odd ; 

.•.I I sin'^ — ^-^ — ^cos2(^ — ^)cos2^d/Adl^ 

= 47r cos 2^' or 0, as i is even or odd ; 
/. cos2<^' 

^ Q 2 , . ,^£?P,(C0S^) ^., 
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in 2 . . ft /v d'/l (cos 0') ^,, 

5.6.7.8 a/ir' ^ 

^ } 

= 2 -^^ 2^ (0±4 + 3X576 -^ 5T6X8 + •-] ' 

ff ence the potential of a spherical shell, of radius c and 
irface density cos 2^', will be 

md 

''^«°^2^' (iX3^ J-«+37i;x6r-+d^P+ ••••) ' 

at an internal and external point respectively. 

19. We will now explain the application of Spherical 
Harmonics to the determination of the potential of a homo- 
geneous solid, nearly spherical in form. The following 
investigation is taken from the Mdcanique Celeste, Liv. iii. 
Chap. II. ft 

Let r be the radius vector of such a solid, and let 

r = a -f a (ttj Fj 4- ttg Fj + ... + a< F, + ...), 

r being a small quantity, whose square and higher powers 
nay be neglected, a^, a.^,...a,... lines of arbitrary length, and 
:^, y,,... F,... surface harmonics of the order 1, 2,...t... re- 
pectively. 

The volume of the solid will be -^ ira^, 

o 

For it is equal to 

== .. Tra', since I / 1\ dfid<f> = 0, 
«> » J -1 . 

jr all values of i 
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Again, if the centre of gravity of the solid be taken as 
origin, a^ = 0. 

For if z be the distance of the centre of gravity from the 
plane of a?y, 



4 /•*• ri rair 

,^T:a*z = j I I r*fidrdfjLd<l> 



-\-ia\{a^y\+a^Y^'\-... 4- c^Y^ ■{-... )dfidij> 



-lUr; ^ 



Similarly 



4 , 
3^ 



^ = ^a" a . ttj I I (1 - /Lt')i cos <f> Y^ dfidj), 

4 _ ri raw 

-7ra'^ = 4{a'a. a^ I I (l-ft')isin<^Y'jd/A6?<^. 

Now Y^ is an expression of the form 

J/A + J5(l-)iA^4cos<^+ C(l-/i,')isin^, 

and therefore all the expressions x, y, z cannot be equal to 0, 
unless ttj = 0. 

We may therefore, taking the centre of gravity as origiBi 
write 

r = a + a{aj^Y^+ ... -^ aiYi+ ...), 

as the equation of the bounding surface of the solid. 

Now this solid may be considered as made up of a homo- 
geneous sphere, radius a, and of a shell, whose thickness is 

The potential of this shell, at least at points whose leas* 
distance from it is considerable compared with its thickness, 
will be the same as that of a shell whose thickness is aa, and 
density 



''•S^*+-+-:^>H' 
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pp being the density of the solid. Therefore the potential, 
for any external point, distant R from the centre, will be 

4nrp,^ + 4>7rp,'xa [-5- ^ + --+-^i:;:i ^i + -"■) ■ 

The potential at any internal point, distant B from the 
oentre, will be made up of the two portions 

I irpJP + 2irp, (a' - J?) or 2,rp U - ^j 

for the homogeneous sphere, 

for the shell, and will therefore he equal to 
2irp, (a'-g-j + 47rp.aa (^^» ^ + ... +^^ -^, + ....) . 

20. If the solid, instead of being homogeneous, be made 
up of strata of diflFerent densities, the strata being concentric, 
and similar to the bounding surface of the solid, we may 

deduce an expression for its potential as follows. Let - r be 

the radius vector of any stratum, p its density, r having the 
same value as in the last Article, and p being a function 
of c onlj. Then, Sc being the mean thickness of the stratum, 
that is the diflFerence between the values of c for its inner 
and outer surfaces, the potential of the stratum at an ex- 
ternal point will be 



'^ILiWr'^"'') ^^^' 



To obtain the potential of the whole solid at an external 
point we must iuteorrate this expression with respect to c, 
between the limits and a, remembering that p is a func- 
tion of c, 

1— "i 
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Again, the potential of the stratum, above considered, 
at an internal point will be 

"*":ii+l c* 7 ^'^' 

To obtain the potential of the whole solid at an iBtemal 
point we must integrate the expression (1) with respec?t to c 
between the limits and R, and the expression {2J witk 
respect to c between the limits R and a, remembering ii 
l>oth cases that p is a function of c, and add the resolti 
together. 



CP. 



CHAPTER V. 



SPHERICiX HARMONICS OF THE SECOND KIND. 

1. We have already seen (Chap. il. Art. 2) that the 
lifferential equation of which P^ is one solution, being of 
kbe second order, admits of another solution, viz. 

•j p: (1 - /*') • 

Now if ^ between the limits of integration be equal 
to ± 1, or to any roots of the equation P, = (all of which 
roots lie between 1 and — 1), the expression under the 
integral sign becomes infinite between the limits of inte- 
gration. We can therefore only assign an intelligible 
meaning to this integral, by supposing /x to be always be- 
tween 1 and 00 , or between — 1 and — oo . We will adopt 
the former supposition, and if we then put (7= — 1, the 

Of 1 \ 

expression p^n _ «) (^^'^' p» (^»,.lJ ^i^l ^® always posi- 
tive. We may therefore define the expression 



P 



/. 



dfi 



p^H-'-iy 

as the zonal harmonic of the second kind, which we shall 
denote by ^„ or Q^ {fi), when it is necessary to specify the 
variables of which it is a function. 

It will be observed that, if fi be greater than 1, P, is 
always positive. Hence, on the same supposition, Q^ is 
always positive. 

Weseethat Q,= -^ — r^^^^g^ir^* 

Jfi /I — L ^ fl»— X 
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d/t 



n r ^1^ 



= 7^ /i lofir r — 1. 



And, in a similar manner, the values of Q^, Qs«*** niay 
be calculated. 



2. But there is another manner of arriving at these 
functions, which will enable us to express them, when the 
variable is greater than unity, in a converging series, with- 
out the necessity of integration. 

This we shall do in the following manner. 

Let J7= , V being not less, and /i not greater, than 

unity. 

™^„ dU_ 1 dU 1_ 

dv~ if-p,)" dii.~{p-iJiY' 

n *\^^ 1^-1 « t,dU l-A* 
IJn _ ^ ^ _ Zui r_2f^ 2 Vgl-M" 

dv\^ '^' dv]~[y-fir)V-l v-(*)~ (v-fi)*' 

d f,, ,,^dtn_ 1-^' / 2/* 2 Ngl-zty. 

Now, let be expanded in a series of zonal harmonics 

p — /I 

P». P,(ji)...P, 0), so that 
^=7^^ = ^» W ^» + ^.M ^. W + ... + *,(") i'.0*)+- 
Then|^|(l-,i')-^^ = ...-,(,-+l)^.WP.(,*) + ... 
Zjf the definition of P (jt) . 
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.d.l»^{(i-^f}-...4|(i-^)^}p,M^... 

And these two expressions are equal Hence, equating 
le coefficients of P,(ft), 



1 
dif 



'(^-''')-dr}=-*"(*+i)^'«- 



Hence 4>i(y) satisfies the same differential equation as P, 
id Q,. But since U^- when i^ = oo , it follows that ^, (i') =0 
hen i/ = 00 . Hence if>^{v) is some multiple of Q^{v)=AQ^{v) 
ippose. It remains to determine A, 

Now, <t>i{v) may be developed in a series proceeding by 

1 
icending powers of - , as follows. 

We have = -+ ^ + ,,,4-_^4-,,„ 

id also =^oW-Po(a^) + <^iW^iW+..- + ^iMP,(a^) + ... 

Now, by Chap. ii. Art. 17, we see that, if m be any 
teger greater than i, the coefficient of P, in /a"* is 

-». ^x (m — 1+ 2) (m — i + 4)...(m— 1) . « . , ,, 
li + 1) -7 — -^-^ — , , , ^ . — tn — 7 — AM — sn if * be odd, 
^ (m + i + l)(m + t-l)...(m + 4)(m + 2) 

— 1 being always even. 
Hence, writing for m successively i, i + 2, * + 4, ... we get 

,(•;) = {li + 1) |(-27:j:i)-(2i-l),..(i + 2) i?^* 

4.6...(t + l) 1 



(2«+3)(2i+l)...(i + 4) v' 



.1+* 



■^(2 



6.8...(i + 3) 1 1 'f 'VkA 1^ 

i + 5)(2t + 3)...(i + 6| v*"^' i 
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2.4...t 1 



and =(2;+l) 



+1 



if t be even. 



4.6...(l + 2) _1_ 

■*" {2i+8) {2i+ l)...^i + 3) v'** 

6.8...(t + 4) 1^ 

■^ (2t + 5) {2i + 3) ... (I + 6) J/**' 

Now, recurring to the equation 

J,.P,(^);»(^*-1) 
we see that, if. Qi{v) be developed in a series of ascending 
powers of -, the first term will be ^ . ,. — ., . , . , , where G 
is the coeflScient of fjf in the development of P, (/a) ; 

thatis C=(^±|l|±i)-vi!;'rl)if,-beodd. 

^. tJ»0»«« (!'""■ Lj 

and = (i±l)(i+3)>-+ 5),...(2^-J) ^^ . ^ ^^^^ 

/* . 4 . 0«**% 

Hence the first term in the development of Q, {y) is 

i2.4. O...!?- ~~ J.) •/••! IT 

;- -1. — ,„\ — ., ; — ^ . - ■ if t be odd, 



(i + 2) (i + 4). . . (2i - 1) (2i + 1) 

and = .. . ,> ,. — «. — >^ .- ... .^ ; — ^pr II t be even, 
(i+1) (i + 3)...(2t-l) (2i-f-l) 

which is the same as the first term of the development of 

P/H divided by ^^—y. 

Hence A = 2i+ 1, and we have 

_ L^= <?,(„) p,0*) + 5Q,{v) P» + 5Q,iv) P.(,*) + ... 

3. The expression for Q, may be thrown into a more 
convenient form, by introducmg m\,o \Xi^ TLxrccL^tator and de- 
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lominator of the coefficient of each term, the factor neces- 
lary to make the numerator the product of i consecutive 
ntegers. We shall thus make the denominator the product 
tf i consecutive odd integers, and may write 



5.6.7...(t + 4) 1 
5.7.9...(2i+5)i/* 



+ li e? 7^ /c\ • . ^^ ~i+6 H" • • • 



(2A?-^l) (2A? +2)..;(^-^2fc ) 1 
"^ (2iH- 1) (2A; + 3)...(2iH- 2A; + 1) i/*^^-'"^ ••" 

whether i be odd or even. 

4. We shall not enter into a full discussion of the pro- 
perties of Zonal Harmonics of the Second Kind. They will be 
ound'very completely treated by Heine, in his Handhuch der 
'{ugelfunctionen. We will however, as an example, investi- 

dQ^ . 

du 



ute the expression for -j-^ in terms of $,^j, ^,+3... 



Recurring to the equation 

+ (2i + l)5Q.(,;) />» + ... 
e see that 

Now we have seen (Chap. 11. Art. 22) that 

Hence ^^i^ = (2t+ 1) P,(/i) + (2i-3) P^(m) +... 
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^^^ = (2i + 9) P^i,.) + (2i+ 5) P^(^) 

+ (2i+l)P,(/A) + ... 

And therefore the coeflScient of P^Ot*) in the expansion 
. d 1 . 

of J 18 

a/M V" fJL 

(2t+l) {(2i+3) Q^,M + (2t+7) <?,«W + (2i+ll) (2^M +...). I 
Again, 

And ^J_ + ^_i_ = 0. 

Hence, comparing coeflScients of P, {fi), 
^^ = - (2i + 3) Q.„ M - (2i + 7) <2^ {v) 

Hence it follows that 

iiQM _ i±M^ =_ (2i + 3) Q^.M, 

and therefore that 

5. By similar reasoning to that by which the existence of 
Tesseral Harmonics was established, we may prove that there ' 
is a system of functions, which may be called Tesseral Har- 
monics of the Second "Kmd, dLeivs^di ixotCL T('> in the same 
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manner as Q^ is derived from P^, The general type of such 
expressions will be 






r/')(/*)l>»-i) 

and this when multiplied by cos a^ or sin (r<f>, will give an 
expression satisfying the differential equation 

{^^ -'*'^ T^ ?7+K(i+i) (1-/*') -a*) tr=o, 

and which may be called the Tesseral Harmonic of the 
second kind, of the degree i and order cr. 



CHAPTER VI. 



ELLIPSOIDAL AND SPHEROIDAL HARMONICS. 



1. The characteristic property of Spherical Harmonics 
is thus stated by Thomson and Tait (p. 400, Art. 637). 

"A spherical harmonic distribution of density on a spheri- 
cal surface produces a similar and similarly placed spherical 
harmonic distribution of potential over every concentric 
spherical surface through space, external and internal." 

The object of the present chapter is to establish the ex- 
istence of certain functions which possess an analogous pro- 
perty for an ellipsoid. They have been treated of by Lam^, 
in his Lefons sur Us fonctions inverses des transcendantes et 
les fonctions isothermes, and were virtually introduced by 
Green, in his memoir On the Determination of the. Exterior 
and Interior Attractions of Ellipsoids of Variable Densities, 
(Transactions of the Cambridge Philosophical Society, 1835). 
We shall consider them both as functions of the elliptic co- 
ordinates (as Lam^ has done) and also as functions of the 
ordinary rectangular co-ordinates; and after investigating 
some of their more important general properties, shall pro- 
ceed to a more detailed discussion of the forms which thejr 
assume, when the ellipsoid is a surface of revolution. 

2. For this purpose, it will be necessary to transform 
the equation 

into its equivalent, when the elliptic co-ordinates €, v, v are 
taken as independent variables. If a, 6, c be the semiaxes 
of the ellipsoid, the two sets of independent variables are 
connected by the relationa 
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J ^ J- = 1 L _ ^ I -s 1 






a'^ + i;' ' b'+v ' c^-f-i/ 

Thus a' + e, J' + e, c' + e are the squares on the semiaxes 
of the confocal ellipsoid passing through the point a;, y, z. 

c? + u, V-\- u, <? + u, the squares on the semiaxes of the 
confocal hyperboloid of one sheet. 

«* + V , V + v\ (? + u , the squares on the semiaxes of the 
confocal hyperboloid of two sheets. 

Thus, e is positive if the point x, y, z be external to the 
given ellipsoid, negative if it be internal. 

And, if a? be the greatest, c' the least, of the quantities 

€ will lie between — (? and oo , 

^ » it ^ >} ^ f 

V „ „ -a' „ -6*. 

3. Now -j-^ + ^-j + -7-7 = IS the condition that 

taken throughout a certain region of space, should be a mini- 
mum. In the memoir by Green, above referred to, this 
expression is transformed into its equivalent in terms of a 
new system of independent variables, and the methods of the 
Calculus of Variations are then applied to make the resulting 
expression a minimum. We shall adopt a direct mode of 
transformation, as follows : 

Suppose a, ^, y to be three functions of x, y, z, such that 

V^a = 0, v')8 = 0, vV = (1), 

wich also that the three families of surfaces represented by 
the equations a = constant, )8 = constant, 7 = constant, inter- 
•ect each other everywhere at right angles, *\.^, %\3l^ NJwbNk 
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dxdx dy dy dzde * dxdx dydy dadz" * 

d^^^d2d0^d2d0^^ 

dxdx dydy dzdz ^ 

Then 

dV dVdoL dVd^ dVdy 
dx d% dx dfi dx d/^ dx* 

dx'''l[l[dx) "*" dff'Kdx) ■*■ dy' \dx) 

^ d'V d^dy g d'Vdydx g d'V doidp 
dfidydxdx dydzdxdx dadfidxdx 

^dV^ dV(Pff dV^ 
doL da? dfi dj? dy da? ' 

d^V d'V 

-j-T Q^d -T-j being similarly formed, we see that, when th 

three expressions are added together, the terms involvin 

dV dV dV 

-T- , -777, -r- will disappear by the conditions (1), and thos 

involving ^ , ^ , ^ by the conditions (2). Hene 



"^ d^ 



+ 



l(D'-(D'-(S)} 



drf 
4. Now, let 

d^ 



J. Ua» 



K<»* + tXi«'"-v ^U^' A -V^X* ' 
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All these expressions satisfy the conditions (1), for a is 
the potential of a homogeneous ellipsoidal shell, of proper 
density, at an external point, and )3 and 7 possess the same 
analytical properties. 

Again, a is independent of v and v, and is therefore con- 
stant when € is constant. Similarly )8 is constant when v is 
constant, and 7 is constant when v is constant. Hence a, fi, 
7 satisfy the conditions (2). 

Now 



U) ■*■ w "*■ U) 



And ■-r:r" + p^+-2-r~ = l- 

f a-' ' f ^ \d€_ 2x 

with similar expressions for ,- and -r- . Hence, squaring 
and adding, 

But from the equations 

• /3 is a purely imaginary quantity. We may, if T^e please, "write \^- IjS' 
for/}. 
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we deduce 

^ ai* ly* ^' _ (® "" ^) (® — v) (« — V ) 

ft) being any quantity whatever. For this expression is of 
dimensions in &>, e, v, v , it vanishes when © = 6, v, or v, 
and for those values of cd only, it becomes infinite when 
ft) = — a*, — i*, or — c", and for those values of ft) only, and it is 
= 1 when ft) = X . 

From this, multiplying by a' + ft), and then putting 
w = — a', we deduce 

,^ (64-a')(i; + a')(i^'4-a') 

a result which will be useful hereafter. 

Again, differentiating with respect to co, and then putting 

ft) = 6, 

x^ y' z^ ^ (e — v) (e ~ v) 

... (*y H. (*)■ + (*y . 4 (i±,«-)ii+ei'+£!) , 

\idx/ \dyj \dz J (e — w) (e — v ) 

•'• fcj ■*■ (d^j "^ fe j " (e - .;) (e - w') ' 

[v - V) {v - e) {e - v) { ' d'-i^ ^ dpT 

The equation V* V= is thus transformed into 



Dr 



+ 
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(v'-v-) [{(e+a*)(e + 6') (6 + c«)}4|JV 
+ (e - v') [{(v + a") (v + 6') (« + 01* I;] V 

(v-e) [^{(„'+a«)(v'+J«)(t,'+c')}4^]V. 

5. A class of integrals of this equation, presenting a close 
analogy to spherical harmonic functions, may be investigated 
in the following manner. Suppose jB to be a function of e, 
satisfying the equation 

he + a') (e + F) (e + c^]^ ^J^= i'^ + ^> 

m and r being any constants. 

Then, if H and H' be the forms which this function 
assumes when v and v are respectively substituted for €, 
the equation V'F= wiU be satisfied by V=^EHH'. 

6. We will first investigate the form of the function 
denoted by E, on the supposition that j& is a rational integral 
function of e of the degree n, represented by 

We see that 

+ + P.} 

= n ["(« - 1) (e + a') (e + i") (e + c') je""* + (n - 2) ;),6"-» 

(n-2)(n-3) 

^ (.+t^ (fH!^ + («+«•) (.^■a^)^■(.^- J) (.+}•) f .„ ^ ^.^ jj^^^„ 
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Hence writing 
we see that 



8' 



n 



Un - 1) (6» + 3/;6« + 3/36 +/0 {e- + (« - 2) p.« 



.«-« 



+ 1 (e« + 2/,e+/J |e"- + (n-l);,.6«^ + ^^ 



(n-2)(n-3) ^ 

(n-l)(«-2) ^ 
2 ft* 



Hence, equating coefiScientsof like powers of e, we get 

« fm + 2 j = m, 

n [^(n - 1) [{n -^)p,+ S/j + 1 {(n - 1) i>, + 2/,}] = nmp, + r. 

n (n - 1) 
= 12 > Pi + ^^i 



f 3 

or, as they may be more simply written, 



=rp. 
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n 



Hn - 1) \7h - ^j p^ + 3n/| = nmp^ + r. 



» {^^^^^^ (« - i) i-. ^ 3 (» - »y^. + K-S/'} 



n(n — 1) 
= — j-^^mp. + nrp, 



n 



|(^ - ^)fzPr^2 + 2/«^«-l} = ^/^« 



It thus appears that p^ is a rational function of r of the 
first degree, p^ of the second, p^ of the n***, and when the 
letters Pj,pj...p^ have been eliminated, the resulting equa- 
tion for the determination of r will be of the (w + 1)*** degree. 
Each of the letters p^^ P%"-Pn ^^ have one determinate 
value corresponding to each of these values of r; and we 

have seen that wi = w(w4-^). There will therefore be (w + 1 ) 

values of E^ each of which is a rational integral expression 
of the ri^ degree, n being any positive integer. 

7. But there will also be values of JE?, of the n*^ degree, 
of the form 

(,+5.)J(e+c^)i|e--+(n-l)g,6-*+(i^:^^ 
We thus obtain 



de 



(e + a*)i (e + 6')(e + cO(n - 1) je^ + {n-i)qjr* 



fn-2)(n-3) ^^ ^ 

j-g J,€ + ... + q^ 



\ 
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+ (6 + a»)i (6 + 6')} {n - 1) je"^ + (n - 2) ?,€""» 

f« _ 9.^ /■»! _ a^ 



+ 



(n-2)(n-3) 

+ 1^ ?1« + — + ?■.! 

(e+o')* (€+**) (e+c*) (n-1) (n-2) |€"-»+(n-3) q,< 

^ (n-3)(n-4) ^ 

Hence 



3)("-4) >. . 



\h 



/^ 1 ^ L«-« .L ^^ 9^ ^ ^*-8 J- (^ — 2) (n — 3) »^ , , 

+ (e+ a') (e+ 6') (e + c") (n-1) (»-2) |€«^+ (n-3) ?,€""* 

(n-3)(n-4) ' 

H j^^ q,e +...+jn 

c 

= (Twe + r) ] 6""* + (n - 1) 2,€""» 

(w-lHw-2) 
+ J— 2 ?»« +•••+?■-) 

.'. (n-1) (2+n-2J = m, 

(n-l)|2a« + |(6' + c')+|(n-2)j.| 

+ (n-l)(n-2)(a* + J'' + c'' + (n-3)gJ = (n-l)«i5,+n 



r» 
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By a similar process to that applied above, we shall find 
"ftliat T is determined by an equation of the v^ degree, and 

t m=(n— 1) f 71 — -j, and that each of the letters j,, 

i--?«-i ^s ^ rational function of r. Thus, there will be n 
lutions of the form 



There will also be n solutions of a similar form, in which 
•be factors (e + c*)* (e + a')*, (e + a*)* (e + fc*)* are respectively 




volved. Hence, the total number of solutions of the v^ 
^degree will be 4n + 1. 

8. We may now investigate the number of solutions of 

degree n-i-^^ n being any positive integer. These will 

!» of the following forms : three obtained by multipljdng a 

3ational integral function of e of the degree n by (e + a^K 

, (e+fc*)*, (e + c*)^, respectively, and one by multiplying a 
'.-ntional integral function of e of the degree n— 1 by the 
l^oduct 

An exactly similar process to that applied above will 
shew us that there will be n + 1 solutions of each of the 
first three kinds, and n of the fourth. Hence the total number 
of such solutions will be 3(n+l) + n, or 4m + S, that is 



hi) 



+ 1. 



To sum up these results, we may say that the total 
number of solutions of the n^ degree is 4?» + 1, n denoting 
either a positive integer, or a fraction with an odd numerator, 
and denominator 2. 

Similar forms being obtained for H, IT, we may proceed 
to transform the expression EHH' into a function of x^ y, z. 

9. Consider first the case in which 



118 ELLIPSOIDAL AND SPHEROIDAL HARMONICS. 

Write this under the form 

E= (e - o)J (e - o)J ... (e - oj. 
Then -0"= (v — Wj) (v — oj^) ... (v— coj, 

ir' = (v'-o,J (v - ft),) ... (i/' -c»„). 
Hence 
EHH'={e-a>^) (v-oj (u -o)J...(€-o)J (v-a)n) (v -©.). 

Now we have shewn (see Art. 4 of the present Chapter) 
that (e — 6>i) (v — fi)j) (u — (o^ 

= (a>, + a0(a), + 6^(a>, + c0f^^ + ^+-^-ll. 

Each of the factors of EHH' being similarly transformed, 
we see that EHH' is equal to the continued product of all 
expressions of the form 

(a, + a^)(a> + y)(a, + c*)f-^ + ^+-5r?^ iV 

^ '^ '^ ^ \a\-\-to 6+0) (T + o) / 

the several values of a> being the roots of the equation 

« .-I ^ (^ "~ 1) — • ^ 

G)" + n;?,©"^^ + \^ ^ py^ + . . . +l>n = 0. 

As this equation has been already shewn to have (w + 1) 
distinct forms, we obtain {n + 1) distinct solutions of the 
equation V'F=0, each solution being the product of n 
expressions of the form 



a* + CD 6' + o) c* + o) 

That is, there will be w + 1 independent solutions of the 
degree 2n in x, y, z, each involving only even powers of the 
variables. 

10. To complete the investigation of the number of solu- 
tions of the degree 2w, let us next consider the case in which E 



f •+...+?, 



«-l 
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The object here will be to transform the product 

(€ + 6«)*(v+y)i(v+J')*(€+c")*(v + 0*(«^'+0^ 

ftince the other factors will, as already shewn, give rise to the 
product of 71 — 1 expressions of the form 

a+o 6' + o) r + o) 

Now, by comparison of the value of a? given in Art. 4, 
We see that 

(6 + 6*) {v + b") {v' + b') (€+0 (v + c^ (v' + O 

= (fc'-c») (i'-a*) (c'-a') (c*-5')3^;5». 

Hence, we obtain a system of solutions of the form of 
the product of (n — 1) expressions of the form 

a? y» z" ^ 

a*+o) J'+o c* + ft> 

multiplied by yz. Of these there will be w, and an equal 
number of solutions in which zx, xy^ respectively, take the 
place of yz. 

Thus, there will be 4?^ + 1 solutions of the degree 2n in 
the variables of which n + 1 are each the product of n 
expressions of the form 



a* + « 6' + 






n are each the product of (n— 1) such expressions, multiplied 

byy«, 

ft ••• ••• ••• *^» 

Tl ... ••• ••• *^» 

11. We may next proceed to consider the solutions of the 
degree 2n + 1 in the variables x, y, z. 

Consider first the case in which 
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Here the product (e + a*)* (v + a*)* {v + a*)* will, as just 
shewD, give rise to a factor x in the product EHH\ 

Hence we obtain a system of solutions each of which is 
the product of n expressions of the form 



a* + ® 6* + a> c* + ft) 



-1, 



multiplied by x. Of these there will be w + 1, and an equd 
number of solutions in which y, z^ respectively take 
place of the fgu^tor x. 

Lastly, in the case in which 



(n-l)(n-2) 
1.2 



.«-« 



p.« +—+i'« 



we see that in EHW the product 

(e+o*)* (w+a*)* (u'+o»)* (e+ J*)* («+&»)* (u' +&«)4 (e+c*)* 

will give rise to a factor xyz. 

Hence we obtain a system of solutions each of which is 
the product of (n — 1) expressions of the form 

sr* V* «* , 

a' + ft) 6* + ft) c^' + ft) ' 

multiplied by xyz. Of these there will be «. 

Thus there will be 4n + 3 solutions of the degree 2» + 1 
in the variables, of which 

(n + 1) are each the product of n expressions of the form 



a? 



.V* 



a' + ft)^ i^ + ft)"*'^^:^"^"''^*'^^^^^^'^' 

(n + 1) are each the product of n such expressions, multiplied 

byy, 
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n are each the product of (n — 1) such expressions, multi' 
iplied by xyz. 

12. Now an expression of the form G . EHB^y G being 

any arbitrary constant, is an admissible value of the potential 

a? V* «' 
mt any point within the shell hj + '^ + -, = 1. But it is 

not admissible for the space without the shell, since it 
liecofnes infinite at an infinite distance. The factor which 
leoomes infinite is clearly E, and we have therefore to 
enquire whether any form, free from this objection, can be 
Edund for this factor. We shall find that forms exist, bearing 
Ihe same relation to E that zonal harmonics of the second 
kind bear to those of the first. 

Now considering the equation 

nrhich we suppose to be satisfied by putting U= E, we see 
that, since it is of the second order, it must admit of another 

^articular integral. To find this, substitute for U, E Ivde, 

we then have 

[{(e+a«)(e+J»)(e + c»)}i^]cr 

= [{(e + a*) (e + t") (e + ,^}i ^ E . jvde 

•- 1 {(e + **) (e + c*) + (e + c*) {.e + a*) + (e + a») (e + b*)} Ev 
^e + a')(e + i^(e + cO(^.v + E^. 
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Now, since by supposition, the equation for the determi- 
nation of Z7 is satisfied by putting Z7= E, it follows that 

when E\vde is substituted for Uy the terms involving tdc 

will cancel each other, and the equation for the detenmna- 
tion of V will be reduced to 

Irft; idE 1/ 1_ . _L_ . _1_ > A 

whence log v + 21og ^ + log {(e + a*) (€ + 6*) (e + c*)}* 

= log t;^+ 2 log -6:^+ log ahc, ^ 

%\ and E^ being the values of v and E, corresponding to €=(llj 

E} abc 



Hence v = v. 












We may therefore take, as a value of the potential at 
any external point. 

For this obviously vanishes when €=oo. It remains so 
to determine % that this value shall, at the surface of the 
ellipsoid, be equal to the value G.EREHy already assumed 
for an internal point. This gives 

C^v^.E^abc \ r. 

Hence, putting », . E* . abc = F,, we see that to the value 
of the potential 

r^ESH'T ^ 



J E* V'^6 -V a*^ V^ ^ b*\ ^^ Jr <? \\4 ' 
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tOr any internal point, corresponds the value 

for any external point. 

13. We proceed to investigate the law of distribution of 
density of attracting matter over the surface of the ellipsoid, 
corresponding to such a distribution of potential. 

Now, generally, if Sn be the thickness of a shell, p its 
volume density, the difference between the normal compo- 
lents of the attraction of the shell on two particles, situated 
ilose to the shell, on the same normal, one within and the 
)ther without will be 4i7rpBn. This is the attraction of the 
ihell on the outer particle, minus the attraction on the inner 
>article. 

But the normal component of the attraction on the outer 

dV 
>article estimated inwards is — 7- • 

an 

And, if F' denote the potential of the shell on an in- 
ternal particle, the normal component of the attraction on 

dV 
t estimated inwards is — ^— . 

an 

Hence ^wpm = ~j -j— . 

'^ an an 

dV^dVdx dVdy dVdz 
dn dx dn dy dn dz dn' 

dec 

And -T- is the cosine of the inclination of the normal at 
dn 

he point a?, y, « to the axis of a?, and is therefore generally 
squal to e-^ — , e denoting the perpendicular from the 
entre on the tangent plane to the surface 



a'-f€ 6*+€ c*-V€ 
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And we have shewn that 

(a» + €)(a» + i;)(a« + i;') 



(^ 



whence 



2e?af_ 1 
X de a' + e' 



a + € a€ 

da? ^ da? 
.'. -y- = 2e-^. 

Similarly -^ = 2e ^ , -^ = 2e -^ , 

•^ a/i ae an ae 

• ^=26^—— ^^ + ^^^ = 26 — 
' ' dn \da? de dy de dz dej de ' 

dV dV 

Similarly -y— = 26 — j— . 

an de 

Now r^V,.EEH'r * 

' Jo ^MCe + a' 



•• de -^-^^ rf€Jo£-{(6 + a")(€ + i»)(6 + c*)li* 
And F= F,.^^^'/^ ^Ka^ + e)(/+e)(c« + e)}i^ 



r d€_ 



therefore, generally, 

de ® de 

But, when the attracted particle is in the immedia 
neighhoMihoo^ of tlie sxntace, €=^. ^^^^^, >(k^^ <Ks%t lii 
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the value of -~ becomes identical with the value of —y— , 

cLe ae 

id we have 

de de ' E, abc' 

> denoting the value which E assumes, when e = 0. 

TTTJ' 1 
Hence, 4,rp8„ = 2e F. -^ ^ . 

But hn, being the thickness of the shell, is proportional to 

and we may therefore write f" = ^ » Sa being the thick- 

iss of the shell at the extremity of the greatest axis ; 

_V, a 1 HE' 
•*' '^"27rSa a6c E, ' 

Qd this is proportional to the value of F corresponding to 
ny specified value of e, since EE' is the only variable 
wtor in either. 

Hence functions of the kind which we are now considering 
ossess a property analogous to that of Spherical Harmonics 
uoted at the beginning of this Chapter. On account of 
his property, we propose to call them Ellipsoidal Harmonics, 
Jid shall distinguish them, when necessary, into surface and 
olid harmonics, in the same manner as spherical harmonics 
re distinguished. They are commonly known as Lamp's 
^^unctions, having been fully discussed by him in his Legons. 
\q equivalent expressions in terms of a?, y, z have been con- 
idered by Green in his Memoir mentioned at the beginning 
f this chapter, and for this reason Professor Cayley in his 
Memoir on Prepotentials," read before the Royal Society 
n June 10, 1875, calls them " Greenians." 

We may observe that the factor 

-1-^ JL 

47r ha abc 

1 11 
equal to . , ^ , aud therefore also to i kt qi -t — v^ - 
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Hence, it is equal to 



4i7r ' 

-s- {hcBa + caBb + abSc) 

or to 1 . 

volume of shell * 

and the potential at any internal point 

/•oo y 

= i volume of shell x JE& . p | r, 

and the potential at any external point 

= i volume of shell x EE^ . p ( — .; 

where for p must be substituted its value in terms of t/andu. 

14. We will next prove that if F^, F, be two different 
ellipsoidal harmonics, dS an element of the surface of the 

ellipsoid, lleV^V^d8=^0, the integration being extended all 

over the surface. 

We have generally 

And throughout the space comprised within the limits of 
integration, V* V^ = 0, V V^ = 0. Hence 



//«(''.f-''.f')^«-»- 



Now it has been shewn already that Fj, F, are each of 
the form EHIT, where JF is a function of € only, H the same 
function of v, R' of v. We may therefore write 

»".=/.(e)/.W/.(v'). 
and aiiailar]/ T^«=/»(«)/»(f)/»(f'). 
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fence v^=Vv£^ 

Lence r, ^ - v, k. ^^ ^^^ , 

* de ~ ^^/ (e) * 
I ow, all over the surface, e = 0. Hence 

/•' fo) f (0) 
[ence, unless *^J-77yf - J (J = 0, which cannot happen 

IS the functions denoted by yj and j^ are identical*, or 
diflfer by a numerical factor, we must have 

jjeV,V,d8=0. 

Tow e is proportional to the thickness of the shell at 
)oint. Calling this thickness 8e, we have therefore 



jjSeV,V,d8=0. 



;e, adding together the results obtained by integrating 
ssively over a continuous series of such surfaces, we get 



jjjV,V,dxdydz = 0; 



V^ now representing solid ellipsoidal harmonics, and the 
ration extending throughout the whole space comprised 
n the ellipsoid. 

This may be shewn more rigoroasly by integrating throngh the 
bounded by two confocal ellipsoids, defined by the values X and /a of e. 
en get, as in the text, 

JJ ^ *(/»(/*) /iW /a(X) M^)S 

be factor within { } cannot yanish for aU yalnes of X and fi, nnless the 
3n9 devoted by fi and /^ be identical, or only differ by a numerical 
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15. It will be well to transform the expression 



//■ 



eV.V^dS 



to its equivalent, in terms of v, v. 

For this purpose we observe that if ds, da be elements of 
the two lines of curvature through any point of the eUipsoid, 
dS^dsd^. 

Now, 
ds" is the value of etc* + dy^ + d!s* when e and v are constant, 
djsl^ ••• ••• ... 6 and u 

and a? - (^ + «') ('^ + «') C'^' + «') . 

therefore if e and v' do not vary, 

^dx dv 

.'. dx^-^-^ dv. 

1 ?/ 1 2? 

Similarly dy= ^ ,/^ dv. dz = ^ -r. dv ; 

.', da^^dx^ + dy^ + dz^ == -r \-r-i t*+ ,lJ^ x« + / / J ^^'- 



Again, diflferentiating with respect to <o the expression 

a-a yi ^ 

obtained for —. h -r^ V -r 1, we get 



i^ , y* g' ^ (v - o)) (v - ft)) 

(a' + a))'"^(6'+ft))*"*"(c' + ft))* Ca' + ft))(*' + a>j ^c^ + o)) 

(v-Q>)(6-ft)) (6-ft))(v--ft))(u -ft)) . 

■*■ ('a'+fi)) (6*^ + 0)) (c^-\-toVQ^"-V^"(6"H- ft)) (c" + ft)) ■*■'••' 
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therefore, putting ai = v, 



aj" 



+ 



y' 



s + 



_^ (v — v ) (e — t;) 



(«• + y)" ^ (i" + v)' ^ (c* + «)* ~ (a» + «) (6* + u) (c» + u) ' 



/. &» = 4 



(y'-v)(e-v) 



di^. 



A similar expression holding for ds* we get 

Jl (v'-v)'(6-v)(f-v) , , , „ 

16 (o'+u) (i'+u) (<;•+«) (a*+u') (b'+v) {(f+v) ' 



Again, 



a> 



+ 



s + 



_ {e-v')(,e-v) 



(o«+6)(y + 6)(c'+e)' 
fvriting e for o> in the expression above ; 

(a*+e)(b'+e)((f+e)iv'-vy 



7z di/dv^. 



•• •^'"'^ 16 (aVi;) (6'+i;) (c»+i;) (a'+v) (6*+i/') (c'+i;') 

It has been shewn that, integrating all over the surface, 
tbe limits of v are — c* and — fi", those of v , — b* and — a*. 

Hence, F^, F,, denoting two different ellipsoidal har- 
^jnonics 



l-v^J'i^Ka^ 



V^rAv-v)dudv' 



[(ff^-v) (y+i;) (c'+v) (o'+i;') (6'+i;0 (c"+i;')}* 



= 0. 



The value of the expression jjjV^dxdjfdz, or its equiva- 



lent 



•^I'-X' 



F" (v'-v)dvdv' 



•» {(a»+w) (y+v) (<;•+«) (o*+ w) (6*+ v') i<f+ v')]* ' 

in any particular case, is most convenientlj obtained bv 
expressing F afi a fiuiction of x, y, z. 

w.B. ^ 
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16. Before proceeding further with the discussion of ellip- 
soidal harmonics in general, we will consider the special case I 
in which the ellipsoid is one of revolution. We must enquire | 
what modification this will introduce in the quantities wbdch 
we have denoted by a, ^8, 7, viz. 



dslr 



^=r ^ 

^f''' d± 

and in the differential equation 

We will first suppose the axis of revolution to be the 
greatest axis of the ellipsoid, which is equivalent to supposing 
A' = c'. To transform a and y, put 0* + -^ = ^, o' + 6=ij', 
a' -\- v = ct>* ; we then obtain 

To transform /8, we must proceed as follows. 

Put i/r = - c' cos' TD- - 6' sin" isr, v = -c'cos*^-6»sin"(^, 
we then get generally 

6» + ^=(6«-c')cosV, c' + ^=(c"-6*)sin*«r; 

cf-^ = 2 (c' — 6*) cos-srsintD-dtsr; 
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:enco, . |..-1(,._„.+ 5,|, 

Iso, € = rf — a\ u = CO* - a*, v = - 6', and our differential 
ion becomes 

(a«-6«_a,»){(,*-a' + &')|JV 
+ (^'-a' + J*) |(a,«-a» + J»)^r F 

-(a,«-,')(a'-5')^=0, 

(»» - a* + J*) |(,» - a* + J*) j-l* r 

_(a»-6')(,»_a,')^^=0. 

bis equation may be satisfied in the following ways. 

irst, in a manner altogether independent of <^, by sup- 
^ F to be the product of a function of rj and the same 
on of 6), this function, which we will for the present 
e hjf(rf) or/(a)), being determined by the equation 



^ {(»'-«• +^')^}/w '-/(»)• 



V-^ 
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if r 

Seoondly, by supposing ^ts & constant multiple of T 
= — 0^V, suppose. 

Our equation may then be written 

-<^(f^-V)[{m*-(^+V)-{v'-a* + b')} F=0, 

which may be satisfied by sajqposiog the factor of V in 
|N;ikdent of ^ to be of the form F{ii) JP{«»), where 

' y- a? + J*) ^lV(i,) - a* (a? - 6*) ^(i,) = m (i|»-a'+}') JJ 

I 

The £Bu;tor inTolving ^ will be of the form 

^ cos o*^ + ^ ffln o*^ 

Now, retaining to the equation 

we see that, supposing the index of the highest power ( 
involved in.f(j^) to be t, we must have m=^i(t-^ 1). 

Now, it will be observed that tf may have any ^ 
however great, but that m\ which is equal to a* + v, n 
lie between o' — 6' and 0. Hence, putting «*= (a*- P 
where /a* must lie between and 1, we get 
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Hence this equation is satisfied by /{(«'— 6')* m} = CP,, 
' being a constant ; and supposing (5 = 1 we obtain the 
)IIowing series of values for/ (o)), 

i = 0, /(a)) = l, 

t = l, y'(a))= 1, 

3o)'-(a'-y) 

1 = 3, /{c») = ^^ — i — -, 

Exactly similar expressions may be obtained for/(i?), and 
lese, when the attraction of ellipsoids is considered, will 
jply to all points within the ellipsoid. But they will be 
ladmissible for external points, since rj is susceptible of in- 
3finite increase. 

The form of integral to be adopted in this case will be 
3tained by taking the other solution of the differential 
][uation for the determination of f{v)t i*®* ^^^ zonal har- 

lonic of the second kind, which is of the form Q, \ j[ , 

here 

d0 






I 



(^-a«+J") 



Or, putting 17* = (a" — h*) i/", ^ = (a* — J') \*, we may write 

d\ 



q,m=p,m/ 



17. We may now consider what is the raeanin? of the 
iiantities denoted by rj and o). They are the values of ^ 
hich satisfy the equation 
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and are therefore the semi-axes of revolution of the 
confocal with the given ellipeoid, which pass through 
point X, y, z. One of these surfaces is an ellipsoid, 
its semi-axis is tf. The other is an hyperboloid of two sk 
whose semi-axis is <o. 

Now, if 5 be the eccentric angle of the point «, jf, 
measured from the axis of revolution, we shall have 

a? = 17* cos* 0. 

But also, since 17*, «•, are the two values of ^ wl 
satisfy the equation of the surface. 

Hence ©• = (a" - 6*) cos" 0, 

and we have already put 

whence the quantity which we have already denoted by 
is found to be the cosine of the eccentric angle of the 
Xf y, z considered with reference to the ellipsoid co: 
with the given one, passing through the point a?, y, z. 
have thus a method of completely representing the poteni 
of an ellipsoid of revolution for any distribution of dei 
symmetrical about its axis, by means of the axis of 
lution of the confocal ellipsoid passing through the 
at which the potential is required, and the eccentric 
of the point with reference to the confocal ellipsoid. For 
any such distribution can be expressed, precisely as inthft 
case of a sphere, by a series of zonal harmonic functions of 
the eccentric angle. 

18. When the distribution is not symmetrical, we must 
have recourse to thefonn of solution which involves the fiM5t(» 
^1 cos <r<^ + jB sin <r<^. It will be seen that, supposing ¥ to 
represent a function of the degree », and putting w = t (i+l), 
the equation which determines F{m) is of exactly the same 
form as that for a te^^ral ^i^hericsJ harmonia For -F(i;), if 
the point be wVWiVsi ^ii^ ^S^^^^j^, ^^ ^^\. ^^^^i ^bsa^ form, 



ET.LTPSOIDAL AND SPHEROIDAL HARMONICS. 135 

• without, representing the tesseral spherical harmonic by 

, or T/*^^ (i/), we adopt the form 



ua — 



(a - 6«)*J 



TfK^lj 



19. It maybe interesting to trace the connexion of sphe- 
cal harmonics with the functions just considered. This may 
e effected by putting i' = a*. We see then that 17 will become 
jual to the radius of the concentric sphere passing through 
le point, and 17* — a' 4- J' will become equal to rf. Hence 
16 equation for the determination of/ (17) will become 

hich is satisfied by putting /{t}) =rf, or rf^^^^K The former 
)lution is adapted to the case of an internal, the latter to 
lat of an external point. 

With regard to /(o)), it will be seen that the confocal 
yperboloid becomes a cone, and therefore a> becomes inde- 

uitely small. But u, which is equal to 1, remains 

uite, being in fact equal to - or cos 0. Hence /(/i) becomes 
16 zonal spherical harmonic. 

Again, the tesseral equations, for the determination of 
'(17), F{a)), become 



( 



v*^)'F{v)'=i(i+l)v*F{v), 



hich are satisfied by F{rD =17* or rf^^^^K 
And, wilting for cw*, {a^-b^) fi\ we have, putting F(ci)) =xif^)> 



{ 



hich gives xW- T^^ (f). 



I 
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20. We will next consider the case in which the axis of 
revolution is the least axis of the ellipsoid, which is equi- 
valent to supposing a* = 6'. To transform a and ft put 
c' + '^ = ^*, c* + 6 = r}\ c* + V = ft)', we thus obtain 

J^ a -c +(f (fl*-(^)i n 

o^oT d6 2 ^ -I a> 

To transform 7, we must proceed as follows : 

Put .'^ = — a'sin'-sr — J'cos'«r, v = — a'sin*^ — J*co8'^, 
we then get, generally, 

a* + i^= (a" - V) cosV, V + 1^= - (a* - 6*) sin'w, 

c'+'^=c*— a'sin*(D-ycos*^ rf'^=— 2 (a*— 6^ sin«r cos«T(2cr. 

Hence 

^ = 2r ^= .^ ?* ifa« = 6'. 

J 4> (a* sinV 4- b* cos*«r — c*)* (a* — c*) « 

Hence, |- = _ ^ (a«_c' + ,») |, 

also, e = 17* — c', 



8 8 

V = 0> — C, 

' 8 



and our differential equation becomes 



V 
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We will first consider how this equation may be satisfied 
•y values of V independent of <j>. 

We may then suppose Fto be the product of a function 
f rj» and the same function of ©, this function, which we will 
appose to be of the degree i, being determined by the 
quation 

^{(a'-C + ,')|}/(,) = f(i + l)/(,), 

On comparing this with the ordinary diflFerential equa- 
Lon for a zonal harmonic, it will be seen that, on account 
f a* being greater than c", the signs of the several terms in 
lie series for /(rf) will be all the same, instead of being 
Iternately positive and negative. We shall thu3 have 

* = 1. Av) = " 



i=%f{v) = 



3i?''+o'-c' 
2(a*-c') ' 



-' ^ 2(a'-c')* 

*-*' -^W 8 (a' - &)* ' 

ad generally 

We will denote the general value of/ (97) by |>, ■! — ^ — --vj- , 
r, writing 17 = (a" - (?)^v, by p^ (y). 
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For external points, we must adopt for/ (17) a function 
which we will represent by qA , , x f > or y,(i/), which will 
be equal to 

It is clear that /(«) may be expressed in exactly the 
same way. But it will be remembered that 17* and «' are 
the two values of y which satisfy the equation 

Hence 97, as before, is the semi-axis of revolution of the 
confocal ellipsoid passing through the point {x, y, z). But 
if CO* = — (a* — c') a", an essentially negative quantity, since 
a* is greater than c*. Hence ©* is essentially negative. Now, 
if ^ be the eccentric angle of the point (x, y, z) measured 
from the axis of revolution, we have a* = 17* co^O. Hence 

,7V = -(a" -017" cos* ^, 

and therefore ©' = — (a* — c*) cos" 

= — (a* — c*) ^J?y suppose. 

Hence the equation for the determination of f{<o) assumes 
the form 

the ordinary equation for a zonal splierical harmonic. Hence 
we may write 

/(o,)=P,W, 

II being the cosine of the eccentric angle of the point a?, y, z, 
considered with reference to the confocal ellipsoid passing 
through it. 
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21. We have thus discussed the form of the potential, 
corresponding to a distribution of attracting matter, sym- 
metrical about the axis. When the distribution is not 
symmetrical, but involves <^ in the form A cos a<f>'\- B sin o-^, 
we replace, as before, P^ {fi) by T^^ (fi), and p^ (fi) by a 
function t^^'^^iv) determined by the equation 



t^'Kp)=(i+vy^Mv), 



1) 



and J, M by t}'H.)[^^^^ 

22. As an application of these formulae, consider the fol- 
lowing question. 

Attracting matter is distributed over the shell whose 

surface is represented by the equation -7 + ^—u — = 1> so 

that its volume density at any point is P, (^), fi being the 
conine of the eccentric angle, measured from the axis of 
revolution ; required to determine the potential at any 
point, external or internal. 

The potential at any internal point will be of the form 

(7P»P.W (1), 

and at an external point, of the form 

CPMQM (2), 

where (a* — J*)* v = the semi-axis of the figure of the con- 
focal ellipsoid of revolution passing through the point (/a, v). 

Now the expressions (1) and (2) must be equal at the 

surface of the ellipsoid, where v = —- — — -. . 

(a* - b'p 

Hence 
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Bat generally 
Hence 



^'V-6«)4~^U(a--6';4j/- 



(a«-6«)* 



■■°'''{i^^-'^^-y^}f. 



dX 



a J>/\\\*r\* 



(O^-ftS)* 



P,(\)I'(V-1) 



We may therefore, putting C = AP^ \—^ ^[ , write 

and we thus express the potentials as follows : 

^Pi (/*) -P« (y) Qi I 2^2 i t ** ^ internal point, 

-4-P« (m) Qi W ^O TT^Ti^jr ^* ^ external point 
Or, substituting for Q^ its value in terms of P,, 

at an internal point, 

at an external point. 

Now, to determine A, we have, Sa being the thickness 
of the shell at the extienxit^ ot Xk^ ^^^ Qlx^^<il\ition, 
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^ 4!irSa.'ri\dri drf J^-a 

4ir Sa o" - i' Uv ~ dv h'-fz^ 

1 



1 a -A. r>/ \ 



^•{(^rr^ijl (a«-j'-^) 



a»-6' 



-1 



~ TZ. wsl -^ "i 0*)' 



4ir 6'&i 
Hence, if p = P, (/»), we obtain 

^ = 47r— =4^J8J. 
a 

And we thus obtain 



.(\)|»(X»-1) 



r, = «» P, w i., {j^.} P, w f ^=^^^ 

If the shell be represented by the equation 
it may be shewn in a similar mannei that 'we ihall bivt^^ 
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V, = 4iwadaP, (fi) q,{v) p, | _^ ^ ■ | . 

23. We may apply this result to the discussion of the 
following problem. 

If the potential of a shell inthe form of an ellipsoid of 
revolution about the greatest pma^ be inversely proportional 
to the distance from one focus, find the potential at any 
internal point, and the density. 

If the potential at P be inversely proportional to the 
distance from one focus 8, and H be the other focus, we have, 

HP+8P = 2v, HP-8P=2(o, 

Hence if M be the mass of the shell, V^ the potential at 
any external point, 

M 



K- 



M 



M 



(a" - 6')4 



S(2i + l)P.0*)<2.(./). 



Now, by what has just been seen, the internal potential, 
corresponding to P, (ji) ^, {y), is 

Hence, if V. be the potenVAai. a.\. «k3 \}CiVetMiSL\Ri\sA, 
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I And the volume density corresponding to P, {fi) Q^ [v) is 



^bSbP. 






Hence the density corresponding to the present distri- 
bation is 

« - ^ ■* /«--.-.x Pi(l*) 

If "F, had varied inversely as JJP, we should have had 

F=^ 

and our results would have been obtained from the forgoing 
by changing the sign of »^ and therefore of fi. 

24. Kow, by adding these results together, we obtain 
the distributions of density, and internal potential, corre- 
sponding to 



J| — i# f; + » ff-^df* 



or, in geometrical language, 

r-j^ Jtf;_^fiP+FP 

= If multiplied by the axis of revolution of the owiSyaJ 
eUipeoid, and divided by the square on the conjugate A^^i- 
diameter. We may express tfau by saying that tbi^ potesr&£;hi 
atanypaiiii aa the dbpaoid is inreisely ytif^n^iyaosiL v> ^^*i 
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square on the conjugate semi-diameter, or directly as the 
square on the perpendicular on the tangent plane. 

Corresponding to this, "we shall have, writing 2k for i; 
since only even values of i will be retained, 

'^ 47r (a" - 6')HS6 p « I 

h being 0, or any positive integer. 

Again, subtracting these results we get 



Jf M ,. 2© 



= -^-^-772. 



CO iy — ft) 



= M multiplied by the distance from the equatoreal plane, 
and divided by the square on the conjugate semi- diameter. 

This gives, writing 2k + 1 for % 



Q: 



M-m 



25. In attempting to discuss the problem analogous to 
this for an ellipsoid of revolution about its least axis, we see 
that since its foci are imaginary, the first problem would re- 
present no real distribution. But if we suppose the exteroal 
potential to be the sum or difference of two expressions, each 
inversely proportional to the distance from one focus, we 
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obtain a real distribution of potential — in the first case 
inversely proportional, to the square on the conjugate 
semi-diameter, in the latter varying as the quotient of the 
distance from the equatoreal plane by the square ou the 
conjugate semi-diameter. 

It will be found, by a process exactly similar to that just 
adopted, that the distributions of internal potential, and 
density, respectively corresponding to these will be ; 



In the first case 



2M 



'Ai^^ 



i being 0, or any positive integer. 
In the second case 



9M" ^*" if^" /.»\5i 



p 2^__2(4fc + 3)^^-^ 



47r(a'-c')iaSa"''"' ' "' „ f <> ]* 

i being 0, or any positive integer. 

26. We may now resume the consideration of the ellip- 
soid with three unequal axes, and may shew how, when the 
potential at every point of the surface of an ellipsoidal shell 
18 known, the functions which we are considering may be 
employed to determine its value at any internal or external 
point. We will begiji by considering some special casea^ 
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by whicli the general principles of the method may be made 
more intelligible. 

27. First, suppose that the potential at every point of 
the surface of the ellipsoid is proportional to a = — ^ suppose, 

Ia this case, since x when substituted for V, satisfies the 

equation v*F= 0, we see that V^- will also be the potential 

at any internal point. But this value will not be admissible 
at external points, since x becomes infinite at an infinite 
distance. 

Now, transforming to elliptic co-ordinates 

And the expression 

V^ f (e+a«)(v+a')(v'+a'') U r" ^2^; 

a \ (a'-b')ia'-c') j j, (^+a»){(Vrta')(^+i*)(t+c')li 



J a 



satisfies, as has already been seen, the equation y*F=0,i9 

equal to V^- at the surface of the ellipsoid, and vanishes 

at an infinite distance. This is therefore the value of the 
potential at any external point. It may of course be written 

r^ r d^ 



*/ 



28. Next, suppose that the potential at every point of 
the surface is proporlioTisl \jo ya=V*^*^^ suppose. In this 
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ime, as in the last, we see that, since yz when sahstitated 
■r V, satisfies the eqnaticm \^V= 0, the potential at any 

^lemal point will he V^^; while, suhstituting for y, z their 

168 in terms of elliptic co-ordinates we obtain {or the 
Ltial at any external point 

if djr 



< 



29. We will next consider the case in which the po- 

iiitial, at eveiy point of the sar&ce, yari^ as ^^ F^^ 

impose. This case materially differs from the two just con- 
id^redy for since af does not, when suhstitated for V, satisfy 
lie equation v'F= 0, tiie potential at internal points caonot 
b general be proportional to ^. We have therefore first to 
avestigate a function of x^ y, «, or of €, v, v which shall 
Mtisfy the equation v*F=0, shall not become infinite within 
klie surface of the ellipsoid, and shall be equal to ^t:' on its 



Kow we know that, generally 

- ((^ + ») (&•+») (<? + »> = (€- ») (v-a>) (v - an). 
And, if 0^, 0^ be the two values of i» which satisfy the 
lunation 

'ipe see that 

and ^{e-0^(y-0;)(v-0^'=O. 

And, by properly detetmiiuiig the ooeflicients ^,, ^j, -4,, 
It is possible to make 

» ^ when 6'c'*' + c'o'y* + o'i V - <^^ - 0. 



X4St JEXLIPSOIDAL AND SPHEROIDAL HABMONICSI, 

Hence, the expression (2) when -4^, A^, -4, are properly 
determined will satisfy all the necessary conditions for an 
internal potential, and will therefore be the potential for 
every internal point. 

Now, we have in general 

-'{a' + 0^(b'^Y0^){<^ + 0^)^{e-0,)(y-0^W-0;^ 

and, over the surface 

fc'c V + cV/ + o'6 V - o»6V =;= 0. 

Hence, ^ being any quantity whatever, we have, all over 
the surface, 

+ 1^;^] (^-^0 C"-^.) (v'-^J -^(^-^J (^-^J, 
and therefore, putting Sr = — a', 

Hence, the right-hand member of this equation possesses 
all the necessary properties of an internal potential. It 
satisfies the general differential equation of the second order, 
does not become infinite within the shell, and is proportional 
to x^ all over the surface. 

We observe, by equation (1), that 
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^denticaUy, and therefore, writing — o* for to, 

Hence, over the surface of the shell, 

and we therefore have, for the internal potential. 

This is not admissible for external points, as it becomes 
infinite. at an infinite distance. , We must therefore substi- 
tute for the factor 6 — ^^ 

(e^0) r. __^ *t_. i— ^ 

with a similar substitution for € — ^^, thus giving, for the 
external potential, 

. r d>lr 



df 

mi 



_^ r d^ 
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The distribution of density oyer the surface, correspond, 
ing to this distribution of potential, may be investigated by 
means of the formula 

'^ 27rda\d€ de J^^o' 
or its equivalent in Art. 13 of this Chapter. We thus find that 

,, 1 g F. r (v-g,)(v--g.) 
'^ ivda Saic I M^^ {0^ - 0^) {a* + 0,) 



(v-gj (v-g,) . f d^lr 

30. The investigation just given, of the potential at an 
external point of a distribution of matter giving rise to a 
potential proportional to a? all over the surface, has an in- 
teresting practical application. For the Earth may be re- 
garded as an ellipsoid of equilibrium (not necessarily with 
two of its axes equal) under the action of the mutual gravi- 
tation of its parts and of the centtrifugal force. If, then, 
V denote the potential of the Earth at any point on or with- 
out its surface, and ft the angular velocity of the Earth's 
rotation, we have, as the equation of its surface, regarded as 
a surface of equal pressure, 

/. F -f- ^ fl* (a? + ^) = a constant, IT suppose. 

Hence, if a, h, c denote the semi-axes of the Earth, we 
iavOj for the determination oi V^^k^fcUamn^ c<mditions: 
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55^ + ^ + 1? = ° W* ■ 

F— at an iufioite distance .(2), 

F= n - 5 n» (a? + y») when 



2 



a* 6* c ^ ' 



The term 11 will, as we know, give rise to an external 
potential represented by 

1 1 

The two terms —^ flV, ^^Cl^y^, will give rise to terms 

which may be deduced from the value of F", just given by 

1 1 

successively writing for V^ — ^HV, and — ^fl*6^ and (in 

the latter case) putting 6* for a' throughout. We thus get 

F= In- 1 ii« («•+&«)! f ^- i 






(«-5.) («- 



^r ^± fly <»' , ^M 
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31. Any rational integral function V of x, y, ^, which 
satisfies tli6 e<^uati6n v'F= 0, can be expressed in a series 
of Ellipsoidal Harmonics of the degrees 0, 1^ 2...i in x, y, z, 
Tor if K be of the degree i, the number of terms in V will 
be (* + l)(*+2)(* + 3) ^ jj^^ ^^ condition jfV^O is 

equivalent to the condition that a certain function of x, y, z 
of the degree i— 2, vanishes identically^ and this gives rise 

to ^^ ^-^ conditions. Hence the number of inde- 

o 

pendent constants in 'F is 

(t + l)(i + 2)(i + 3) (i-i)t(/+l) 
6 " 6 ' 

or (i + 1)'. And the number of ellipsoidal harmonics of the 

1 3 t 

degrees 0, 1, 2...i in 4:, y, z or of the degrees 0, 3, 1, ^...^ 

in e, V, V , is, as shewn in Arts. 6 to 10 of this Chapter, 

1 + 3 + 5+.. . + 2i + l^ 

or (i + 1/. Hence all the necessary conditions can be satis- 
fied. 

32. Again, suppose that attracting matter is distributed 
over the surface of an ellipsoidal shell according to a law of 
density expressed by any rational integral function of the 
co-ordinates. Let the dimensions of the highest term in this 
expression be i, then by multiplying every term, except those 
of the dimensions i and z — 1 by a suitable power of 

we shall express the density by the sum of two rational inte- 
gral functions of a;, y, z of the degrees % z^l, respectively. 
The number of terms in these will be 
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Anil any elJipsoidal surface harmonic of the degree %,i — 3... 
in X, y, z, may, by suitably introducing the factor 

be expressed ^s a homogeneous function of x, y, z of the 
degree i; also any such harmonica of the degree i—\, »— 3... 
in X, y, z may be similarly expressed aa a homogeneous 
function of x, y, z of the degree t — 1. And the total number 
of these expressions will, as just shewn, be (i + 1)', hence by 
assigning to them suitable coefficients, any distribution of ' 
density according to a rational integral function of m, y, b 
may be expressed by a series of surface ellipsoidal harmonics, 
and the potential at any internal or external point by the j 
corresponding aeries of solid ellipsoidal harmonics. 

33. Since any function of the co-ordinates of a point o 
the surface of a sphere may be expressed by means of a series 
of surface spherical harmonics, we may anticipate that any 
function of the elliptic co-ordinates u, v' may be expressed by 
a series of surface ellipsoidal harmonics. No general proof, 
however, appears yet to have been given of this proposition. 
But, assuming such a development to be possible at all, it 
may he shewn, by the aid of the proposition proved in 
Art. 1.5 of this Chapter, that it is possible in only one way, 
in exactly the same way as the corresponding proposition 
for a spherical surface is proved in Chap. IV. Art. 11. 

Tlie development may then be effected as follows. De- ] 
noting the several surface harmonics of the degree i in x, y, g, 

or I in u, V, by the symbols V,i'\ V<^\...V,^'-<*^\ and by j 
J^(u, v) the expression to he developed, assume 

F (y, v") = C„ F"„ -t- £7 «)r,(" -I- C«> r,<« -|- C,(« r.W + ...' 
+ C, t" F."> + ... + C,W VM + . . . 

Then multiplying by eF/'* and integrating all over the J 
surface, we have 

jeF (v, v'} yw dS = C.w (e (F,"')' dS. 
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The values of jeF {v, v") Vy^ dS, and of Je ( V^^y dS must 

be ascertained by introducing the rectangular co-ordinates 
X, y, z, or in any other WAy which may be suitable for the 
particular case. The coefficients denoted by C' are thus 
determined, and the deyelopment effected* 
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1. Prove that (sin 0)« = A p^_l|p, + ^i.^. 

Why cannot (sin Gf be expanded in a finite series of spherical 
harmonics % 

111 '^ +SU1 ^ 

2. Provethatl + 2-P,+ 3-P. + ^Pa+— =log r 



sin ;r 

2 



3. Establish the equations 



«i>. = (2n-l)^^.-(n-l)P._.. 

4. If /A = cos tf, prove that 

.^ li + w / fl\* 

!?♦*• i I / 1\* ' f ~i«»« 1 



P.(^) = l-i(i + l)8in'- +... + (- ir^l^^=^(8m'D+.. 



and also that 

^*w=(-iy+(-ir^(*+i)cos»|+,.. 

(1 ^) h-^ \ 2/ 

5. Prove that, if a be greater than c, and i any odd 
integer greater than m^ 



6. Prove that j (^*) V = ».(* + 1). 
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7. Prove that, when m = * 1, -ri; = —* "Wr — • 

^ '^ ' Hit!* i — m 2 \ m 



8. Prove that 






1, Pj ft.* p, 

p p p 

••• *•• ••• ••• 

P P P 



is a nmnerical multiple of 



9. Prove the following equation, giving any Laplace's co- 
efficient in terms of the preceding one : 



^«+i = P^. + ^JJ^-^P+C', 



where Op = fifi+Jl -/x' ^1 -/a" (sOs (« -xd') and (7 is zero if n he 
even, and 

10. If i, y, k be three positive integers whose sum is even, 
prove that 



i 



^ l^S.^O'+Jfe-t-l) 1.3...(A;4-t-i-l) 1.3...(t+i-^-l) 
2.4. ..0* + ^-*) 2.4...(A; + i-j*) 2.4...(i+i-A;) 

2.4...(t+i + ;^) 1 

1.3 ... (i + j + ^-1) i+J+^ + 1' 

Hence deduce the expansion of F^Pj in a series of zonal 
harmonics. 

11. Express a:^i/ + y^ + yz + y + z as a sum of spherical 
harmonics* 

12. Find all the independent symmetrical complete harmonies 
of the third degree and of the fifth negative degree. 

13.' Matter is distributed in an indefinitely thin stratum over 
the surface of a sphere whose radius is unity, in such a manner 
that the quantity of matter laid on an elemfent (SS) of the surface 
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where x, y, z are rectangular co-ordmatea of the element iS re- 
fefi**! to tlie centre as origin, ami a, b, c, f, g, h are constants, 
Pind the value of the potential at any point, whether internal 
or external. 

14. If the radius of a sphere be r, and itE law of density ba 
p = OB + 6^ + cs, where the origin ia at the centre, prove that its 

liotential at aa external point (f, 1;, {) is y--™ («£+ 6ij + c^) where 

R is the distance of (^, i], ^) fi-om the origin. 

\a. I*t a spherical portion of an infinite quiescent liquid bo 
separated from the liquid round it by an infinitely thin flexible 
membrane, and let this membrane he suddenly Hot in. motionj 
every part of it in the direction of the radius and with velocity 
equal to S„ a harmonic function of position on the surface. Find 
the velocity produced at any external or internal point of the 
liquid. State the corresponding proposition in the theory of 
Attraction. 

16. Two circular rings of fine wire, whose masses are MkoA J 
Jr, and rudii a and a', are ])laced with their centres at distances 
, b', from the origin. The lines joining the origin with the 
Jnti'ea are perpendicular to the planes of the rings, and are in- 
elined to one another at an augle 0. Shew that the potential of 
one ring on the other is 



I 



Ind ff^ and Q^ are tiie xatne functions of 6' and a' and of cos 6 

la d sin 6 respectively, and c is the greater of the two quantities 
Va' + 6' and Jd' + b'', 

17. A uniform circular wire, of radius a, charged with 
electricity of line-density e, surrounds an uninsulated concentric 
fpherical conductor of radius c; prove that the electrical density 
at any point of the surface of the couductoi" is 



-iO- 



•i-:^--^i;^- 



tlie pde of the plane of the wire being ihe pole of the haxmonios. 
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18. Of two spherical conductors,, one entirely surrounds the 
other. The inner has a given potential, the outer is at the 
potential zero. The distance between their centres being so 
small that its square may be neglected, shew how to find the 
potential at any point between the spheres. 

19. If the equation of the bounding surface of a homo- 
geneous spheroid of elHpticity c be of the form 



r=a(l-|cP.). 



prove that the potential at any external point will be 

M C-A 



r» 



P,> 



where C and A are the equatoreal and polar moments of inertia 
of the body. 

Hence prove that V will have the same value if the spheroid 
be heterogeneous, the surfaces of equal density differing from 
spheres by a harmonic of the second order. 

20. The equation 2? = a(l+a2/) is that of the bounding 
surface of a homogeneous body, density unity, differing slightly 
in form and magnitude from a sphere of radius a; a is a 
small quantity the powers of which above the second may be 
neglected; and y is a function of two co-ordinate angles, such 
that 

where F^, T ,,,Zq, Z^ .,, are Laplace's 'functions* Prove that 
the potential of the body's attraction on an external particle, 
the distance of which from the origin of co-ordinates is r, is 
given by the equation 
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21. If Jf be the miass of a unifomi hemisplierical shell of 
radius c, prove that its potential, at any point distant r from 
the centre, will be 



3 f^ 3.5 r^ \ 

■*'2.4.6^•c*"2,4.6.8^^c•^••7' 



2r^2 \2^V 2.4 V 



, 3 c^ 3.5 c' \ 

'*'2.4.6^'r«"2.4.6.8 ^V ■*" -7' 

according as r is less or greater than c; the vertex of the hemi- 
sphere being at the point at which fi = 1. 

22. A solid is bounded by the plane of xy^ and extends to 
infinity in all directions on the positive side of that plane. 
Every point within the circle a* + y* = a*, « = Ois maintained at 
the uniform temperature unity, and every point of the plane xy 
without this circle at the uniform temperature 0. Prove that, 
when the temperature of the solid has become permanent, its 
value at a point distant r from the origin, and the line joining 
which to the origin is inclined at an angle 6 to the axis of z will 
be 

r \ r' 1 3 r* 

P - P - 4.-P - P 1. + 

1.3...(2^-l) r^» 

"^^^^ 2.4...2i -^«*+i a«+» + - 

if r < a, and 

1 a^ 1^ a* 1.3...(2t^l) a« 

2^r-"2.4 'r*"*"**'"^ ^ 2.4...2i ^«*+ii:^* + - 

if r > a. 

23. Prove that the potential of a circular ring of radius <?, 
whose density at any point is cos m<^, c\\f being the distance of the 
point measured along the ring from some fixed point, is 



+ ... 
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2. cos m<f> (sin e)'^^-^±-^ ^ -^, 

■^2.4.6...(2to + 2) d/*-" »^ 

1.3.5...(2A!-1) <rf.^^ e'** 1 
*2.4.6...2(«» + A rf/t- j^Htt+i+'-j- 

where r is greater than c« If r be less than e, r and c must be 
interchanged. 

24. A solid 13 bounded by two confocal ellipsoidal surfaces, and 
its density at any point P varies as the square on the perpendicular 
from the centre on the tangent plane to the confocal ellipsoid 
passing through P. Prove that the resultant attractioii.' of such 
a solid on any point external to it or formiDg a part of its mass 
is in the direction of the normal to the confocal ellipsoid passing 
through that point, and that the solid exercises no attra^^on on a 
point within its inner surface^ 



^^. 
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